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FOREWORD 


This  report  covers  the  progress  made  during  the  first  year  of  a 
three -year  program  conducted  by  the  Department  of  Aerospace  Engi¬ 
neering,  The  University  of  Michigan,  Ann  Arbor,  Michigan,  under 
Contract  F08635-71-C-0083  with  the  Air  Force  Armament  Laboratory, 
Eglin  Air  Force  Base,  Florida.  Major  Philip  J.  Unrein  and  Lieutenant 
Richard  E.  Swanson  (DLIF}  served  as  project  monitors  for  the  Arma¬ 
ment  Laboratory.  This  effort  was  begun  on  20  January  1971  and  was 
completed  on  19  January  1972. 

This  research  was  directed  by  Professor  J.  A.  Nlcholls.  The 
contractor  report  number  is  UM  004880 -2-T. 


This  technical  report  has  been  reviewed  and  is  approved. 


ABSTRACT 


This  report  covers  the  progress  made  In  the  first  year  of  a  three - 
year  research  program.  For  convenience,  the  work  Is  divided  Into  two 
phases:  Theoretical  Analysis  and  Experimental  Research. 

The  theoretical  aspects  of  cylindrical  strong  blast  waves  and  Chapman- 
Jouguet  (C-J)  detonations  are  treated  In  detail.  Comparisons  are  made 
between  the  two,  and  a  critical  radius  Is  discussed  which  divides  the  cloud 
into  an  inner  blast  wave  zone  and  an  outer  detonation  zone.  Expressions 
for  the  evaluation  of  ground  and  dynamic  impulse  are  presented.  Other 
aspects,  treated  In  less  detail,  include  propagation  of  a  C-J  detonation 
through  a  cloud  of  fuel  drops,  the  initiation  problem,  and  deviations  from 
the  ideal  cylindrical  blast  model.  Also,  a  modified  computer  program 
for  calculating  detonation  velocity  of  complex  hydrocarbon  fuels  is  dis¬ 
cussed. 

The  design  and  operation  of  an  experimental  facility  to  study  the 
propagation  of  two-phase  detonations  in  a  cylindrical  segment  of  a  cloud 
is  described.  The  generation  of  strong  shockB  and/or  detonation  is 
achieved  by  explosion  of  a  condensed  explosive  (Detasheet) .  The  results 
of  a  number  of  experiments  on  the  propagation  of  strong  cylindrical 
shocks  are  presented  and  compared  with  theory.  The  agreement  is  rea¬ 
sonable  but  some  questions  remain.  Experiments  conducted  with  kerosene 
drops  in  air  showed  that  at  small  radius  the  cylindrical  wave  decayed  as  a 
shock  wave,  but  beyond  a  critical  radius  the  wave  apparently  becomes  a 
constant  velocity  two -phase  Chapman -Jouguet  detonation.  Lower  initiation 
energies  show  evidence  of  a  non-sustalned  detonation.  The  experimentally 
determined  critical  radius  agreed  quite  well  with  theoretical  predictions. 
The  results  lend  encouragement  to  the  prediction  of  threshold  energy  levels 
required  for  detonation  initiation  in  various  geometries. 
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SECTION  I 
INTRODUCTION 


I 


The  aim  of  this  research  program  is  to  gain  a  better  understanding  of 
some  facets  of  the  unconfined  explosion  of  a  distributed  liquid  fuel -gaseous 
oxidizer  mixture.  In  the  general  case  of  such  a  two -phase  cloud,  the  mean 
drop  size  and  fuel -oxidizer  ratio  are  variable  throughout  the  cloud.  Under 
these  conditions  the  ignition  criterion  and  maintenance  of  detonative  com  - 
bustlon  are  much  in  question.  Certainly  it  is  possible  that  detonation  may 
be  established  and  then  quenched  in  other  parts  of  the  cloud  where  the  drop 
size  is  too  large  or  the  mixture  is  too  lean  or  rich.  Another  factor  of 
importance  is  the  finite  size  of  the  cloud  and  the  proximity  to  the  ground 
or  solid  walls.  The  unconfined  boundaries  allow  the  detonation  wave  to  be 
relieved.  That  is,  the  fact  that  the  acoustic  impedance  of  the  surrounding 
gaseous  phase  is  relatively  low  allows  the  pressure  behind  the  detonation 
to  be  reduced  by  lateral  expansion  of  the  gases.  This  relief  leads  to  weak¬ 
ening,  and  possibly  quenching,  of  the  wave.  The  mean  drop  size  will  be 
Important  in  this  sense. 


An  analytical  model  which  is  capable  of  predicting  the  time -dependent 
characteristics  of  a  fuel -air  explosion  is  completely  lacking.  This  is  in 
marked  contrast  to  the  point  explosion  problem  where  one  can  calculate 
the  pressures,  temperatures,  velocities,  etc. ,  for  any  distance  from  the 
explosion  and  for  any  Instant  of  time.  Also,  one  can  immediately  assess 
the  Influence  of  the  density,  heat  of  combustion,  and  properties  of  the 
burned  gases  of  the  explosive.  A  similar  level  of  understanding  is  desir¬ 
able  for  the  fuel-air  explosion,  although  the  problem  is  appreciably  more 
complicated.  Accordingly,  one  phase  of  this  research  program  is  directed 


to  gaining  such  an  understanding.  A  discussion  of  the  approach  used  and 


The  second  general  area  of  study  in  this  program  stems  from  the  fact 
that  controlled  experiments  on  the  cylindrical  or  spherical  propagation  of 
two -phase  detonations  have  never  been  conducted.  The  liquid  fuel  drop 
size  and  the  distribution  of  drop  size  and  fuel -air  ratio  throughout  the  cloud 
are  of  importance  to  fuel -air  explosions,  and  their  influence  would  have  to 
be  estimated  in  the  theoretical  evaluation  just  described.  In  order  to  be 
more  realistic,  an  experimental  program  has  been  initiated  on  a  scaled- 
down  model  of  the  fuel-air  explosion.  Further,  the  laboratory  model  can 
be  used  to  ascertain  initiation  energy  requirements  for  various  mixtures. 
The  results  to  date  from  these  studies  are  described  in  Section  m, 
Experimental  Research. 


SECTION  n 

THEORETICAL  ANALYSIS 

A.  INTRODUCTION 

In  the  establishment  of  a  fuel-air  explosion  (FAE),  liquid  fuel  in  an 
appropriate  container  is  dispersed  into  the  atmosphere  as  a  cloud  of  fine 
droplets  by  the  detonation  of  a  primary  charge.  An  appropriately  placed 
and  timed  secondary  explosion  is  then  used  to  detonate  this  cloud  of  fuel 
droplets. 

The  shape  of  the  fuel  cloud  which  is  generated  is  generally  quite  com¬ 
plex  and  depends  upon  the  design  of  the  primary  charge  and  the  fuel  canister. 
Often  the  cloud  is  doughnut -shaped  with  a  hole  in  the  center  generated  by  the 
munition  wake.  Information  about  particle  size  distribution  and  about  the 
distribution  of  particles  through  the  cloud  is  limited;  however,  it  is  clear 
that  the  cloud  is  far  from  having  a  uniform  distribution  of  monodisperse 
droplets.  The  secondary  blast  may  be  initiated  at  the  center  of  the  cloud,  or 
several  detonators  located  within  the  periphery  of  the  cloud  may  be  used. 

The  detonation  of  the  cloud  may  start  at  the  center  and  move  outward,  but 
in  the  case  of  the  doughnut -shaped  clouds  the  detonation  may  also  propagate 
in  an  azimuthal  direction  around  the  cloud  periphery  and,  of  course,  there 
may  be  more  than  one  center  of  detonation. 

Generation  of  ground  and  dynamic  Impulse  and  of  a  high  temperature 
environment  are  key  objectives  of  the  FAE.  A  goal  of  FAE  research  is 
clearly  to  accomplish  these  objectives  as  effectively  as  possible  with  a 
given  amount  of  tael.  Hence,  it  is  desirable  to  know  the  detailed  mechanism 
of  the  dissemination  process,  a  subject  beyond  the  scope  of  the  present 
research.  Once  a  tael  cloud  is  established,  it  is  important  to  know  the 
strength  of  secondary  blast  required  to  detonate  the  fuel  cloud.  The  man¬ 
ner  in  which  the  detonation  then  propagates  through  the  cloud  will  deter¬ 
mine  the  effectiveness  of  the  FAX. 


In  this  report  the  initiation  and  propagation  of  the  detonation  through  the 
fuel  cloud  will  be  considered.  It  is  unlikely  that  a  detailed  analysis  of  the 
complex  FAE  as  described  above  will  be  possible  so  that  a  complete  under¬ 
standing  must  rely  heavily  on  experiment.  However,  much  can  be  learned 
from  the  analytical  and  experimental  study  of  relatively  simple  models  of  the 
actual  FAE,  and  this  is  the  approach  followed  here. 

As  a  first  step,  an  Idealized  FAE  with  cylindrical  symmetry  and  a  uni¬ 
form  monodisperse  droplet  distribution  will  be  considered.  This  model 
provides  a  starting  point  for  the  consideration  of  other  FAE  phenomena.  At 
very  short  times  the  idealized  cylindrical  FAE  behaves  like  a  cylindrical 
blast  wave,  while  at  long  times  the  behavior  approaches  that  of  a  cylindrical 
Chapman -Jouguet  wave.  These  two  limiting  solutions  are  discussed  in  detail. 

Expressions  for  the  ground  impulse  and  the  dynamic  impulse  for 
Chapman -Jouguet  detonations  and  blast  waves  with  plane,  cylindrical,  and 
spherical  symmetry  have  been  derived.  By  using  the  similarity  solutions 
the  equations  for  Impulse  assume  a  particularly  simple  form  which  make  it 
relatively  easy  to  compare  the  impulse  due  to  blast  waves  and  detonations  of 
various  geometries.  Universal  functions  have  also  been  obtained  which  can 
be  used  to  compute  the  total  Impulse  at  a  fixed  distance  from  the  center  of 
explosion  for  planar  and  spherical  C-J  waves,  and  several  sample  calcula¬ 
tions  have  been  carried  out. 

As  indicated  above,  a  critical  question  is  whether  a  secondary  blast 
will  initiate  a  detonation  in  the  fuel  cloud.  By  reviewing  the  theory  of  blast 
initiation  of  detonations  in  pre -mixed  gaseous  fuels,  it  has  been  possible 
to,  at  least  qualitatively,  establish  some  criteria  for  the  blast  initiation  of 
fuel  droplet  clouds. 

The  idealized  analysis  forms  the  basis  for  studying  important  FAE 
phenomena  which  can  be  treated  as  deviations  from  the  Idealized  model. 

One  Important  effect,  which  is  omitted  In  the  idealized  analysis,  is 


associated  with  the  finite  size  of  the  fuel  cloud.  Thus,  the  high  pressure 
behind  the  C-J  detonation  causes  the  gas  behind  the  wave  to  expand  into  the 
inert  region  surrounding  the  cloud.  The  ground  pressure  therefore  starts 
dropping  to  atmospheric  at  some  finite  distance  behind  the  wave.  The 
effect  of  this  side  relief  is  to  reduce  the  impulse  generated  by  the  FAE. 

In  this  report  the  effects  of  side  relief  are  described  in  the  light  of  current 
work  on  this  subject. 

The  droplet  size  distribution  and  the  distribution  of  the  fuel  withir.  the 
cloud  may  also  have  important  effects  on  the  FAE.  Some  modifications  of 
the  ideal  model  which  will  make  it  possible  to  examine  the  influence  of  a 
non-uniform  distribution  of  fuel  are  considered. 

B.  THE  CYLINDRICAL  BLAST  MODEL 
1.  Description  of  the  Model 

An  Idealized  model  consisting  of  a  cylindrical  wave  propagating  through 
a  monodisperse  cloud  of  uniformly  distributed  fuel  droplets  provides  an 
excellent  starting  point  for  the  theoretical  study  of  the  FAE  process.  The 
secondary  blast  is  initiated  by  the  instantaneous  release  of  energy  Eq  per 
unit  length  of  the  axis  of  the  cylindrical  wave.  This  idealized  model  is 
represented  in  Figure  1.  The  idealized  model  is  equivalent  to  a  cylindrical 
FAE  confined  between  two  non-yielding  parallel  surfaces  (for  n  ‘gligible 
side  relief)  as  also  shown  in  Figure  1.  A  segmented  test  secti  on  is  cur¬ 
rently  being  used  to  simulate  this  idealized  model  of  the  FAE  n  the 
laboratory^.  The  simplified  model  described  above  is  amemble  to  both 
analysis  and  laboratory  experimentation  and  can  also  form  tl  e  basis  for 
the  consideration  of  complications  Introduced  by  the  departure  of  the  actual 

FAE  from  the  ideal.  The  propagation  of  a  blast  wave  into  a  detonating  gas 

(2)  (31 

has  recently  been  discussed  by  Chernyl  et  al  ' ,  Korobeinikov  ,  Bishlmov 
et  al^\  and  Lee  et  al^.  From  these  analyses  it  can  be  concluded  that  at 
short  times  after  the  secondary  explosion  a  strong  blast  wave  will  propa¬ 
gate  into  the  fuel  cloud,  and  at  first  the  energy  released  by  combustion 
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Confining 

Surface 


Symmetry 


Cylindrical  Blast  Wave  Equivalent  Confined  Wave 

Figure  1.  Idealized  FAE  Model 

will  have  a  negligible  effect  upon  the  propagation  of  the  wave.  Provided 
a  detonation  can  be  Initiated  this  blast  gradually  decays  to  a  cylindrical 
Chapman-Jouguet  (C-J)  wave.  It  Is  useful  to  define  a  critical  blast  wave 
radius  i\  such  that  the  fuel  energy  contained  with  r*  Is  equal  to  the  blast 
wave  energy  Eq  [  Korobeinikov*8*  ] .  Thus  r*  Is  given  by  the  expression 

r,  *  (vEQ/ofip^v  U) 

with  v  I' ,  2,  and  3  for  plane,  cylindrical,  and  spherical  geometry.  The 
geometric  factor  oy  Is  given  by 

ov  *  2{v  -  1)  w  +  (v  -  2)(v  -  3) 

Q  la  the  combustion  energy  released  In  the  fuel  cloud  per  unit  mass  of  mix 
ture,  while  pt  Is  the  Initial  density  of  the  fuel  cloud.  Thus,  when  rg  <  r, 
the  Idealised  FAE  behaves  like  a  strong  cylindrical  blast  wave,  and  with 


r  »  r*  the  limiting  FAE  behavior  approaches  that  of  a  cylindrical 
s 

C-J  wave,  where  r  is  the  shock  radius. 

s 

Similarity  solutions  of  the  equation  of  motion  for  the  two  limiting  cases 

of  fuel  air  explosion  are  available  and  are  described  in  detail  by  Sedov. 

No  simple  solutions  exist  in  the  range  r  ~  0(r,)  where  a  transition  from 

blast  wave  to  C-J  wave  behavior  takes  place.  Numerical  methods  or 

(71 

semi -empirical  methods  similar  to  that  proposed  by  Bach  et  al  must  be 
used  to  describe  the  flow  in  this  transitional  regime  which  is  also  intimately 
related  to  the  initiation  problem.  If  r*  «  r  ,  where  r  is  the  radius  of  the 
fuel  cloud,  the  FAE  will  behave  dominantly  as  a  cylindrical  C-J  wave,  and 
this  is  the  regime  in  which  practical  FAE's  will  operate.  Consequently, 
the  blast  wave  and  C-J  similarity  solutions  have,  as  a  first  step,  been 
adapted  to  the  study  of  FAE  behavior  prior  to  undertaking  the  extensive 
numerical  work  needed  to  study  behavior  in  the  transitional  regime. 

A  simple  idealized  model  of  the  FAE  thus  consists  of  the  blast  wave 

solution  for  r  <  r*  and  the  C-J  solution  for  r  >  r*.  Although  this  model 

will  not  be  accurate  near  the  center  of  the  secondary  blast  and  in  the 

(8) 

transitional  regions  near  r  ~  0(r  J ,  one -dimensional  tests  '  and  some 
of  the  measurements  in  the  segmented  test  section  as  reported  below  indi¬ 
cate  that  this  simplified  model  comes  remarkably  close  to  actually  des¬ 
cribing  what  happens  in  a  FAE.  The  strong  blast  wave  and  C-J  wave 
similarity  solutions  thus  play  a  central  role  in  the  study  of  FAE  phenomena 
and  hence  are  described  in  some  detail  below  following  essentially  the 
treatment  of  Sedov!®) 

2.  Similarity  Solution  for  a  Strong  Blast  Wave 

Both  the  blast  wave  and  C-J  solutions  are  based  on  very  general  dimen¬ 
sional  considerations.  In  the  Euler ian  approach  the  dependent  variables 
of  Interest  are  the  fluid  velocity  v,  the  density  p,  and  the  pressure  p.  These 
quantities  will  depend  on  a  linear  space  coordinate  or  radius  r  and  upon  the 


time  t  as  well  as  other  dimensional  constants  which  enter  the  problem 
through  boundary  and  initial  conditions. 

In  the  case  of  a  blast  wave  the  characteristic  independent  parameters 
are 


' 1 


t 


and  y 


where  p^  =  initial  gas  density 
p^  -  initial  gas  pressure 

Eq  =  energy  release,  energy  release  per  unit  length,  or  energy 
release  per  unit  area  depending  on  the  geometry 

y  =  ratio  of  specific  heats. 

From  the  Buckingham  Pi  theorem  it  then  follows  that  the  dimensionless 
velocity  V,  pressure  P,  and  density  R  defined  by 

v  =  (r/t)  V  ;  p  =  PjR  ;  p  =  Pj  (rVt2)  P  (2) 


depend  on  y  and  on  the  dimensionless  ratios 


A  *  f «  « 


1  t2/-2 


(3) 


*2  ~  v/v+2 2/v+2  .-2v/i/+2 

P <  t 


'1 


(4) 


For  strong  blast  waves  such  that  p^/p^  »  1  where p,j Is  the  pressure  be¬ 
hind  the  shock,  it  is  reasonable  to  assume  that  -  0  and  tnen  =  0  so  that 
V,  R,  and  P  depend  only  on  the  single  dimensionless  variable  A.  In  this 

.  i 

case  then 

v  *  (r/t)  V(A)  ;  p  «  Pj  R(A)  ;  p  «  Pj  (r2^2)  P(A) 
so  that  the  solution  Is  self  similar. 


(5) 


In  the  self  similar  blast  wave  the  dimensionless  variable  A  will  have 


a  constant  value  on  the  shock  and  from  this  the  conclusion  follows  that 


r 

s 


t2/e+2 


«5) 


j/o  |  /«  n  /e 

i.e.,  r  varies  as  t  ,  t  ,  and  t  '  for  strong  planar,  cylindrical,  and 

D 

spherical  blast  waves.  Remarkably,  this  result  follows  from  purely  dimen¬ 
sional  considerations  without  any  recourse  to  the  conservation  equations 
and  their  solutions. 


The  flow  behind  the  blast  wave  can  only  be  determined  once  V(A) ,  R(A) , 
and  P(A)  are  known,  and  the  conservation  equations  must  be  solved  to 
determine  these  variables.  The  introduction  of  the  variable 


,  _  v  pW 

z*Ym 


(7) 


in  addition  to  V,  R,  and  P  reduces  the  equations  for  the  conservation  of 
mass,  momentum,  and  energy  to  the  following  set  of  ordinary  differential 
equations. 

dz  z[2(V  -  1)  »  v(r  -  1)  VI  (V  -  0) 
m  [V(V  -1XV-S) 


»{(y  - 1)  v(v  -  i)(v  -  S)  *  [«v  - 1)  *y(r  - «]  z} 
(v  -  a)[v(v  -  t)(v  - «)  +  v  1  -  v  ]  z 


(8) 


<KM  .  z  -  (v  -  _ 

M  V{V  -  1)(V  -  «)  +  v(l  -  Vjz 


(9) 


V(V  -  1)(V  -  «)  +  v 


!B 


*  -  (V  - 


(10) 


where  5  =  2/(v  +  2).  Equation  (8)  is  a  single  first  order  equation  relating 
z  and  V,  and  uncoupled  from  the  other  equations.  The  study  of  blast  wave 
solutions  now  revolves  around  an  analysis  of  Equation  (8)  relating  z  and  V. 
Once  z(V)  is  determined  R  and  P  can  be  determined  by  quadrature  or,  as 
will  be  seen  below,  by  using  integral  forms  of  the  conservation  equations. 
A  study  of  the  solution  curves  in  the  z-V  plane  reveals  many  important 
properties  of  self  similar  solutions  of  one -dimensional  unsteady  gas  flows, 
and  is  carried  out  in  great  detail  by  Sedov^  and  Oppenheim  et  al^ 

The  physical  significance  of  the  variable  z  is  not  always  evident. 
Assuming  that  the  fluid  satisfies  a  perfect  gas  equation  of  state,  it  follows 
from  the  definition  of  P(A)  and  R(X)  that 


z 


(11) 


where  a  is  the  speed  of  sound.  Thus  z  may  be  related  to  the  temperature 
or  to  a  Mach  number.  For  t  4  0,  z  -  «as  r  —  0,  l.e. ,  z  becomes  Infinite 
at  the  center  of  symmetry. 

A  study  of  the  z-V  plane  indicates  that  consistent  solutions  of  the  blast 
wave  problem  do  not  exist  without  the  introduction  of  a  shock  discontinuity, 
i.e. ,  it  is  impossible  to  find  a  continuous  trajectory  in  the  z-V  plane  con¬ 
necting  the  center  of  symmetry  with  the  undisturbed  flow  at  r  -  oc.  Thus, 
the  shock  conditions  in  terms  of  z  and  V  must  be  established. 

The  shock  velocity  c  is  readily  evaluated  since  the  shock  position  cor¬ 
responds  to  a  constat  value  of  A  *  A*.  Hence,  from  the  definition  of  A  in 
Equation  (3),  it  follows  that 


(12) 


(13) 


dr  r  _  r 

_ s  _  c  _  j  _s  _  2  s 


v  +  2  t 


For  a  very  strong  shock  wave  propagating  into  a  stationary  fluid  introduc¬ 
tion  of  the  variables  z,  P,  V  and  R  and  the  Equation  (13)  for  the  shock 
velocity  into  the  Rankine -Hugoniot  equations  leads  to  the  following  result 


for  conditions  behind  the  shock: 


V2  =  (y  +  l)(v  + 


R  ~Y  +  1 
R2  "  y  -  1 


p  = - - - 

2  (y  +  l)(v  +  J)2 


8y(y  -  1) 

2'(y  +  l)2(y*J)2 


Ahead  of  the  blast  wave  and  hence  -  0  and  for  the  strong  blast 
with  =  0,  it  follows  that  z^  =  0  also.  Thus,  the  undisturbed  flow  corres¬ 
ponds  to  the  origin  of  the  z-V  plane,  while  the  point  immediately 

behind  the  shock  Is  determined  by  Equation  (14).  The  problem  now  is  to 


determine  the  variation  of  z  with  V  between  the  shock  and  the  center  of 


symmetry.  Most  treatments  of  blast  wave  theory  at  this  stage  use  numer¬ 
ical  integration  of  the  similarity  equations  to  determine  the  flow  behind  the 
blast;  however,  Sedov  was  able  to  obtain  the  following  analytical  expres¬ 
sion  for  the  function  z(V)  by  ingenious  use  of  the  Integral  forms  of  mass, 
momentum,  and  energy  conservation. 


(y  -  1)  V*  (v  -  irf5' 


9  1 

{v  +  2)y 


that  is  an  explicit  solution  of  z  as  a  function  of  V  is  obtained  and  the  blast 
wave  problem  is  essentially  solved.  Using  the  solution  (15),  for  z(V), 
Equation  (9)  can  be  Integrated  analytically  for  X  as  a  function  of  V.  The 
dimensionless  density  R  can  be  determined  as  a  function  of  z  and  V  using 


integral  forms  of  the  conservation  of  mass  and  entropy  with  the  result1 


,  ...... 


(16) 


■•'’-’('’•rbj' 

where  C ^  is  now  a  constant  of  integration  whose  value  can  be  established 
using  the  shock  conditions  (14) . 

One  key  difficulty  remains.  Although  dimensional  considerations  are 

2/(  iz-t-2) 

sufficient  to  establish  that  r  ~  t  ,  the  precise  variation  of  shock 

s 

radius  r  (t)  is  given  by 
s 

rs  =  x.<yPl> 1/1+2 ‘2M2  (12) 

however,  the  value  of  the  constant  A*  is  not  known.  To  resolve  this  diffi- 
culty,  Sedovv  ’  replaced  the  energy  release  Eq  with  a  quantity  E  oc  Eq 
but  without  a  predetermined  value.  Then  the  choice  of  A*  becomes  arbi¬ 
trary  and  for  convenience  the  value  A*  =  1.0  is  used.  The  relation  be¬ 
tween  E  and  Eq  now  depends  on  the  choice  of  A*  and  can  be  obtained  by 
noting  that  Eq  equals  the  total  energy  contained  between  the  origin  and 
shock.  For  a  cylindrical  blast,  for  example,  equating  Eq  to  the  total 
blast  energy  leads  to  the  following  relation  between  Eq  and  E: 


with  similar  relations  for  planar  and  spherical  symmetry.  Thus,  although 
the  similarity  analysis  suffices  to  etrtabllsh  how  the  shock  radius  rg  varies 
with  time,  the  complete  solution  of  the  blast  wave  problem  is  needed  to 
determine  the  precise  variation  of  rs(t). 


With  =  1,  quantities  at  the  shock  are  given  by  the  following  relations: 
rg  .  (E /p/M2  t2/v*2 


P2  = 

V 


4  /EX172*11  i-v/v+2 

( V  +  2)(y  +  1)  VPj/ 

8pl  /E'?/2+v  t-2v/v+2 
(i v  +  2) 2  (y  +  1)  'pl' 

8(y  -  1)  /E\2/2^  t-2v/v+2 

(v  +  2)2(y  +  l)2  W 


(18) 


The  similarity  variable  A  now  can  be  seen  to  be  the  ratio  of  the  radius  r 
to  the  shock  radius,  i.e. , 

X  -  rAE/pj)171*2  t2/v*2  =  (r/r8)  (19) 


Sedov's  analytical  solution  has  been  programmed  for  the  IBM  360 
computer  so  that  all  parameters  of  blast  wave  flow  can  be  determined  for 
different  values  of  Eq  and  for  different  geometries  with  minimal  computa¬ 
tion  time.  This  program  has  been  used  to  assess  the  early  stages  of  FAE 
and  to  compare  FAE  and  blast  wave  dissemination  of  explosive  energy. 

To  more  clearly  illustrate  the  features  of  blast  waves,  a  computer -animated 
motion  picture  of  the  blast  wave  solution  has  been  prepared. 

3.  The  Self  Similar  Chapm an  - Jougu et  Solution 

As  mentioned  above,  the  C-J  wave  represents  the  limiting  behavior 
of  the  Idealized  FAE  for  large  times.  The  determination  of  the  self 
similar  solutions  for  this  limiting  flow  following  Sedov^  is  outlined  below. 


As  before,  the  variables  of  interest  are  v,  p,  and  p,  and  the  independent 
variables  are  again  r  and  t;  it  is  also  useful  to  introduce  the  dimensionless 
variables  V,  R  and  P  defined  by  Equation  (2).  In  developing  the  self  similar 
solution  the  wave  is  treated  as  a  discontinuity  with  heat  release  Q  per  unit 
mass,  propagating  into  the  unburned  gas  with  constant  speed  c.  The  charac¬ 
teristic  parameters  of  the  problem  are  pit  p1?  Q,  r,  t,  and  y,  and  since 
the  dimensions  of  Q  can  be  expressed  in  terms  of  p^  and  p^,  there  are 
actually  only  4  independent  dimensional  quantities  in  the  C-J  wave  problem. 
Hence,  if  Pj,  Q  and  t  are  chosen  as  the  primary  variables,  only  the  single 
dimensionless  ir  product 


*1  =  y/Qt 


(20) 


can  be  formed  from  the  characteristic  parameters.  The  propagation  of  the 
C-J  wave  is  thus  self  similar  from  the  outset  and,  unlike  the  blast  wave 
problem,  no  assumptions  need  to  be  made  regarding  the  pressure  p^  in  the 
undisturbed  medium  or  the  strength  of  the  wave.  Instead  of  using  as  the 
similarity  variable,  it  is  more  convenient  to  use  the  variable 

X*  (21) 

Vq  t 


where  ^  Is  an  arbitrary  constant  which  is  discussed  below.  Since  A  will  be 
constant  on  the  detonation  discontinuity  in  the  self  similar  case,  it  is  evi¬ 
dent  from  Equation  (21)  that  rfl  ~  t,  i.e. ,  the  wave  will  propagate  with 
constant  speed  for  plane,  cylindrical,  or  spherical  geometry. 

Once  again,  It  is  convenient  to  use  the  variables  z,  V,  R,  and  P,  and 
then  the  conservation  equations  can  be  reduced  to  the  equations 


14 


relating  z  and  V  and  the  equation 


dfln  A)  z-(V-l)2  , 

**  ~v[(V-l)2-w:] 

The  flow  behind  the  detonation  is  assumed  to  be  isentropic,  and  then  the 
integral  forms  of  mass  and  entropy  conservation  yield  the  relation 


with  Cg  a  constant  of  integration.  When  dimensional  quantities  are  intro¬ 
duced  in  Equation  (24),  this  equation  reduces  to  the  statement  p /pY  =  const, 
i.e. ,  the  entropy  is  constant  in  the  region  inside  the  detonation. 

A  continuous  solution  of  Equation  (22)  extending  from  the  center  of 
symmetry  to  the  undisturbed  flow  at  r  -  oc  Is  not  possible  so  that  the 
jump  conditions  across  the  detonation  must  be  Introduced  in  order  to 
determine  the  variation  of  z  with  V.  The  variation  of  \  with  V  can  be 
determined  from  Equation  (23)  by  quadrature  once  the  variation  of  z  with 
V  Is  known.  Unlike  the  blast  wave  problem  discussed  above,  no  analytical 
solution  of  Equation  (22)  for  z  as  a  function  of  V  Is  known. 

For  propagation  Into  a  stationary  gas  with  =  0,  the  Hugoniot  condi¬ 
tions  across  the  detonation,  when  expressed  in  dimensionless  variables, 
become 


(27) 


z,  =  (l  +  -i)  (t  -  A)(l  +  y.  A) 

2  (y2  +  l)2  '  ’l >  2 

The  parameter  A  is  defined  by  the  relation 


(28) 


A  =  0  corresponds  to  a  C-J  detonation  for  which  the  velocity  downstream  of 
the  detonatlve  discontinuity  is  sonic  with  respect  to  the  wave  front.  A>  0 
corresponds  to  a  strong  or  overdriven  detonation  while  A  <  0  corresponds 
to  weak  detonations  which,  from  thermodynamic  considerations,  cannot  exist. 

As  before  A  =  A*  =  const  at  the  det  onative  discontinuity  in  the  self 
similar  case.  It  is  readily  shown  that  choosing  the  convenient  value  A*  =  1.0 
is  equivalent  to 

A  =  r/rfi  (29) 


as  in  the  case  of  the  blast  wave  solution. 


In  order  to  find  Rg,  V2>  and  z 2  from  Equation  (25)  through  (27),  it  is 
necessary  to  know  z^.  Now 


*i 


rtTr.-i-j 

r 

s 


(30) 


2 

where  Is  the  Mach  number  of  the  detonation.  For  a  C-J  wave  Mp 
can  be  determined  In  terms  of  Q,  a^,  y^,  and  y2  from  Equation  (28)  with 
A  «  0. 


The  key  problem  Is  now  to  obtain  a  number  leal  solution  of  Equation  (22) 
starting  from  the  center  of  symmetry  V  *  0,  z  *  ceand  ending  at  the  point 
V  =  V„  z  -  z9  at  the  C-J  discontinuity.  For  this  purpose  it  is  useful  to 
study  the  behavior  of  the  solution  in  the  z-V  plane  shown  In  Figure  2  below, 
where  a  C-J  solution  trajectory  from  the  center  of  symmetry  D  to  the  undis- 
turbed  gas  at  r  -  oc  (point  0)  Is  shown.  The  parabola  z  =  (1  -  V)  corres¬ 
ponds  to  points  where  the  velocity  is  sonic  with  respect  to  the  discontinuity 


and  hence  the  point  (Vg,  Zg)  at  the  C-J  front  lies  on  this  parabola.  The 
point  A  where  V  =  0,  z  =  1.0  is  a  singularity  of  Equation  (22)  and  can  be 
shown  to  be  a  nodal  point.  Hence,  all  C-J  solutions  pass  through  the  point 


A.  It  is  readily  shown  that 


(dr/dt)A  =  a 


(31) 


Figure  2.  Behavior  of  Solution  to  Equation  (22)  In  z-V  Plane 
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In  the  physical  plane  the  point  A  thus  moves  out  with  the  speed  of  sound  a 
and  r  .(t)  =  at  is  a  characteristic  separating  a  stationary  core  of  fluid  from 
an  isentropic  expansion  behind  the  C-J  front.  Physically  the  flow  will  then 
be  as  shown  in  the  graph  of  v  versus  r  below.  In  both  Figures  2  and  3  the 
arrow  is  in  the  direction  of  increasing  X. 

The  detailed  behavior  of  z(V't  near  the  singularity  A  has  been  deter¬ 
mined  analytically  and  has  been  used  in  the  numerical  integration  of  Equa¬ 
tion  (22).  From  Equations  (22)  and  (23)  is  is  also  possible  to  arrive  at 
some  Interesting  results  regarding  the  velocity  gradient  dv/dr  at  the  det¬ 
onation  front  and  on  the  characteristic  r  =  r^.  Since 

X  =  (r/ct)  ;  v  =  (r/t)  V  =  XcV 

it  follows  that 


r 

Figure  S.  Behavior  of  Solution  to  Equation  (22)  in  v-r  Plane 
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At  the  detonation  (point  Hg  in  (he  z-V  plane)  d(ln  X)/dV  =  0  for  v  =  2,  3  so 

that  (dv/dr)  *  oc  immediately  liehlnri  the  detonation  front  for  cylindi  ical 

(5) 

and  spherical  waves.  As  indicated  by  Lee  et  al  ,  this  raises  the  question 
whether  true  spherical  and  cylindrical  C-J  waves  can  exist.  For  planar 
waves  with  v  =  1,  dv/dr  remains  finite  at  the  detonation.  On  the  character¬ 
istic  r  =  r^  it  can  similarly  be  shown  that  (dv/dr) ^  remains  finite  for  plane 
and  cylindrical  waves  but  (dv/dr)  ^  =  0  for  spherical  waves. 

i 

4.  Comparison  of  the  Blast  Wave  and  C-J  Wave  Solutions 

To  show  the  difference  in  the  behavior  of  the  C-J  and  blast  wave  solu¬ 
tions  we  have  compared  these  in  the  planar  case.  Specifically,  the  critical 
blast  wave  radius  r,  which  in  the  planar  case  is  given  by 

r,  =  EQ/2p1Q  (la) 

was  chosen  to  be  20  feet  and  Q  was  based  on  a  C-J  detonation  propagating 
through  a  stoichiometric  methane  (CH^)  air  mixture.  Such  a  mixture  is 
representative  of  a  typical  hydrocarbon -air  mixture.  Knowing  r„,  and 
Q,  the  value  of  Eq  for  the  blast  wave  solution  is  readily  calculated  from 
Equation  (la) . 

For  a  detonation  through  a  stoichiometric  mixture  of  CH^  and  air  at 
Pj  *  atm,  Tj  «  537°R,  Q  is  readily  calculated  following  Elsen  et  al^  with 
the  result  that 

Q«  1203.9  BTU/lbm 
and  then  from  Equation  (la) 

Bq  «  3405  BTU/ft2  *  2. 649  x  106  ft-lb/ft2 
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which  Is  equivalent  to  1. 68  lbm  of  TNT  per  ft  .  A  constant  value  of 
Y  *  1. 3  was  used  in  the  calculations. 


In  Figures  4,  5,  and  6  the  variation  of  p,  v,  and  T  are  shown  for  the 
blast  wave  as  it  propagates  to  the  critical  radius  with  the  corresponding 
quantities  shown  for  the  C-J  detonations  in  Figures  7,  8,  and  9.  It  is 
evident  that  the  blast  wave  and  C-J  wave  flows  are  quite  different.  The 
presence  of  the  quiescent  core  is  evident  in  the  C-J  curves.  The  blast 
wave  pressure  is  initially  higher  than  the  C-J  value  but  rapidly  drops 
to  a  value  about  equal  to  the  constant  pressure  behind  the  C-J  wave.  The 
temperature  T  rises  rapidly  behind  the  blast  wave,  and  T  —  ac  as  r  -  0 
for  the  blast  solution  which  thus  is  not  valid  near  r  =  0.  The  density, 
which  is  not  shown  in  the  figures,  drops  off  rapidly  behind  the  blast  wave, 
and  in  some  approximate  treatments  it  is  assumed  that  the  fluid  in  the 
region  r  <  rg  is  concentrated  in  a  small  region  behind  the  shock. 

C.  CALCULATION  OF  GROUND  AND  DYNAMIC  IMPULSE 

Both  the  ground  Impulse  and  the  dynamic  impulse  are  important 
parameters  in  assessing  the  effectiveness  of  fuel  air  explosions  as  com¬ 
pared,  for  example,  to  a  blast  wave  with  the  same  energy  release.  Both 
the  total  Impulse  imparted  to  the  entire  region  enclosed  by  the  wave  front 
and  the  Impulse  per  unit  area  at  a  point  which  is  a  fixed  distance  from  the 
center  of  explosion  have  been  considered. 

1.  Calculation  of  Total  Impulse 


The  total  Impulse  imparted  to  the  ground  by  either  a  blast  wave  or  a 
C-J  wave  is  given  by  the  Integral 


<e|sd)  ajnssfljd 


21 


Figure  4.  Blast  Wave  Solution  Variation  of  Pressure 


Figure  5.  Blast  Wave  Solution  Variation  of  Velocity 
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(p  -  p.)  rdr  dt 


(34) 


Kt)  -  2./  / 
0  0 


for  cylindrical  and  spherical  waves,  and  I(t)  Is,  oi  course,  a  function  of 
time.  For  the  self  similar  blast  wave  and  C-J  solutions  the  integrals 
(33)  and  (34)  can  be  reduced  to  a  particularly  simple  form. 

For  strong  blast  waves  p^  «  p  and  hence  can  be  neglected  in  (33) 
and  (34).  For  plane  and  cylindrical  symmetry,  l.e. ,  v  =  1,  2,  the  inner 
integrals  in  (33)  and  (34),  when  expressed  in  terms  of  dimensionless 
variables,  become 


The  integral  on  the  right  side  of  Equation  (35)  will  be  a  function  only  of  y 
for  the  self  similar  blast  wave  solution,  hence  letting 


It  is  readily  seen  that  the  grou  id  impulse  due  to  a  strong  blast  wave  will  be 

:  *"1-2  (31> 

that  is,  the  ground  impulse  varies  linearly  with  the  time.  The  universal 
function  0(?)  must  be  evaluated  numerically  from  the  blast  wave  solution. 

For  spherical  blast  waves  the  situation  Is  somewhat  different  for  while 
the  wave  is  spherical  the  ground  Impulse  Is  still  applied  to  a  circular  area, 
as  In  the  cylindrical  case.  Then  It  is  readily  shown  that 
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where  now  P(A)  corresponds  to  the  spherical  blast  solution.  Not  surpris¬ 
ingly,  I(i)  does  not  Increase  as  rapidly  for  the  spherical  blast  wave  as  for 
cylindrical  and  planar  waves.  It  is  important  to  note  that  in  Equation  (37), 
Eq  is  blast  energy  per  unit  area  or  length.  The  impulse  given  by  Equation 
(38)  is  that  generated  by  a  point  explosion  at  the  ground  surface;  however, 
since  the  spherical  blast  wave  solution  applies  to  an  explosion  in  free  space, 
Eq  in  Equation  (38)  is  twice  the  energy  release  of  the  ground  explosion. 

For  planar  and  cylindrical  C-J  detonations  the  ground  impulse  from 
Equation  (33)  will  be 


The  quantity  within  the  brackets  will  be  constant  for  any  given  fuel  oxidizer 

mixture;  thus,  the  ground  Impulse  from  planar  and  cylindrical  detonations 
2  3 

varies  as  t  and  t  ,  respectively.  The  ground  impulse  produced  by  a  hemi¬ 
spherical  C-J  detonation  can  be  obtained  from  Equation  (39)  with  v  =  2  but 
with  P(A)  evaluated  from  the  solution  for  a  spherical  C-J  detonation. 

Equations  (37)  through  (39)  provide  a  relatively  simple  means  of  com¬ 
paring  the  ground  impulse  from  various  types  of  waves.  The  variation  of 
T(t)  with  wave  radius  also  Is  readily  established  from  Equations  (37)  through 
(39).  Thus,  l(t)  oc  r#^*f  rg3,  and  rg3/^  for  plane,  cylindrical,  and 
hemispherical  waves. 

L  % 

Actually  geometric  considerations  indicate  that  the  comparison  of  a 
cylindrical  FAE  to  a  cylindrical  blast  wave  is  not  an  appropriate  basis. 
Thus,  as  shown  in  Figure  10  below,  if  all  the  combustion  energy  in  a 
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Figure  10.  Fuel  Air  Explosion*,  Geometric  Considerations 


cylindrical  FAE  were  concentrated  at  one  point  and  released  instantaneously, 
the  result  would  be  a  spherical  rather  than  a  cylindrical  blast  wave.  Conse¬ 
quently,  the  appropriate  comparison  with  regard  to  ground  impulse  is  be¬ 
tween  a  spherical  blast  and  a  cylindrical  FAE  and  such  comparisons  can  be 
made  using  the  results  developed  above. 


Because  of  the  induced  velocity  there  will  be  a  dynamic  force  on  any 
body  or  obstacle  located  in  the  wake  of  either  the  blast  wave  or  the  detona¬ 
tion.  The  behavior  of  an  object  in  the  path  of  the  wave  is  shown  qualitatively 
in  Figure  11  below.  As  the  shock  sweeps  over  the  obstacle,  It  is  diffracted, 
and  the  flow  during  this  diffraction  process  will  be  unsteady  and  quite  com¬ 
plex.  After  sufficient  time  has  elapsed,  it  seems  reasonable  to  assume 
that  the  flow  will  be  quasi -steady,  that  is,  at  any  instant  the  flow  will  behave 
like  a  steady  flow  with  a  free  stream  velocity  equal  to  v,  the  velocity  behind 
the  shock  or  blast  wave. 

During  the  quasi-steady  portion  of  the  flow  the  maximum  pressure  on 
the  body  will  equal  the  local  stagnation  pressure  pQ  and  will  occur  at  the 
forward  stangation  point.  There  are  now  several  ways  of  estimating  the 
fluid  dynamic  force  on  a  body  in  the  flow.  The  maximum  force  per  unit 
area  will  be  of  the  order  of  pQ  -  p  where  p  is  the  static  pressure  so  one 
estimate  for  the  impulse  at  a  given  position  might  be 
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(40) 


Another  approach  would  be  to  use 
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where  A^  is  the  frontal  area  of  the  body  and  Cp  is  the  drag  coefficient.  If 
the  flow  is  subsonic,  as  for  instance  lx>hlnd  a  C-J  detonation,  it  is  possible 
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in  the  present 
analysis  this 
part  of  flow  is 
neglected 
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Assume  that  after  shock 
passage  the  flow  past  a 
sphere  behaves  as  a 
steady  flow  with  free 
stream  velocity  u  and 
free  stream  pressure  p. 


Figure  11.  Behavior  of  an  Object  in  the  Wake  of  Blast  or 

Detonation  Wave 
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to  assume  isentropic  flow  in  computing  pQ.  Then  if  the  gas  behind  the  wave 
is  assumed  perfect, 


-iV'v1' 


1  + 


M 


where  y 2  is  the  ratio  of  specific  heats  and  M  is  the  Mach  number.  Generally 
[(r2  -  l)/2]  M2  «  1.  Even  when  M  =  1.0,  [{y2  -  l)/2]  M2  =  0. 1.  Then 


P0  2 

1+~MZ  = 
P  2 


In  this  case  Equation  (40)  becomes 


(42) 


Here  Equation  (42)  is  used  as  a  measure  of  Irf,  the  dynamic  impulse  on  an 
obstacle.  Comparing  Equations  (41)  and  (42),  it  can  be  seen  that 

CD  A.  I.  =  Impulse  on  body  with 
drag  coefficient  Cj) 
and  frontal  area  A^. 
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It  should  be  emphasized  that  the  effects  of  the  shock  diffraction  on 
the  Impulse  applied  to  any  body  in  the  flow  Is  neglected  in  using  Equation 
(42)  as  a  measure  of  the  dynamic  Impulse.  The  correctness  of  this 
assumption  should  be  investigated. 

The  impulse  as  written  in  Equation  (4$  represents  the  dynamic 
impulse  felt  by  a  body  at  some  fixed  distance  from  the  explosion  center. 
The  detailed  evaluation  of  this  quantity  shall  be  considered  below.  First, 
however,  measure  of  the  overall  Impulse  induced  by  the  wave  is  needed. 
For  planar  and  cylindrical  waves  the  total  dynamic  Impulse  induced  by  the 


wave  in  a  region  of  unit  height  shall  lie  used  so  that  the  total  dynamic 
Impulse  Td  becomes 


For  a  strong  blast  wave  Equation  (43)  can  be  reduced  to 


which  is  very  similar  to  the  result  for  the  total  ground  impulse.  In  the  case 
of  the  strong  blast  wave,  a  combination  of  Equation  (43)  and  (35)  and  the 
definition  of  a(y)  in  Equation  (17)  leads  to  the  result 


(45) 


so  that  an  expression  representing  what  Is,  in  essence,  a  combined  dynamic 
and  ground  impulse  is  directly  proportional  to  the  blast  energy  and  the  time, 
in  the  case  of  planar  and  cylindrical  waves. 


In  spherical  waves  it  is  more  difficult  to  define  a  total  dynamic  impulse 
which  can  be  related  to  dynamic  effects  in  the  region  enclosed  by  the  wave. 
In  particular,  it  becomes  difficult  to  compare  total  ground  and  total  dynamic 
Impulse  since  the  ground  impulse  is  applied  to  a  circular  area  just  as  for 
cylindrical  waves.  If  pv  /2  is  summed  over  the  volume  of  the  sphere,  the 
total  dynamic  impulse  is  that  given  by  Equation  (44)  with  v  =  3,  or  one -half 
that  value  for  a  hemisphere  bounded  by  the  ground.  The  simple  relation 
[Equation  (45)]  cannot  be  applied  in  the  spherical  case. 
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For  C-J  detonations  the  total  dynamic  impulse  will  be  given  by 
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2.  The  Impulse  at  a  Fixed  Point 

The  ground  and  dynamic  impulse  per  unit  area  at  a  point  a  fixed  distance 
from  the  center  of  explosion  wili  now  be  considered.  In  particular,  the 
ground  impulse  due  to  a  plane  C-J  detonation  travelling  down  a  closed  tube 
and  the  dynamic  impulse  due  to  a  spherical  C-J  detonation  have  been  com¬ 
pleted.  By  making  suitable  approximations,  universal  functions  have  been 
derived  which  can  be  used  to  compute  ground  and  dynamic  impulse,  in  gen¬ 
eral,  for  these  two  configurations. 

In  the  case  of  a  detonation  propagating  down  a  closed  end  tube,  the  key 
boundary  condition  is  that  the  velocity  of  the  fluid  must  be  zero  at  the  end  of 
the  tube.  Consequently,  the  C-J  detonation  is  followed  by  an  isentropic  ex¬ 
pansion  in  which  the  velocity  of  the  fluid  drops  to  zero  as  shown,  for  ex¬ 
ample,  in  Figures  7,  8,  and  9.  The  problem  of  the  detonation  travelling 
down  a  closed  tube  has  been  treated  by  Sedov^  and  Taylor^  ^  among  others. 

In  this  discussion  the  C-J  wave  is  treated  as  an  exothermic  discontin¬ 
uity.  Immediately  behind  the  C-J  wave 


p  *  pressure  =  p^j 

and  the  velocity  of  the  fluid  is  sonic  with  respect  to  the  wave  front,  1.  e. , 

D  ‘  UCJ  *  aCJ 

where  a^j  Is  the  speed  of  sound  in  the  burned  gases  immediately  behind  the 
front,  and  D  is  here  used  for  the  propagation  velocity.  As  before,  the  gas 
behind  the  detonation  Is  assumed  to  behave  as  a  perfect  gas. 


Within  the  isentropic  expansion  liehind  the  wave  the  pressure  is  given  hy 
.  (U) 

the  expression 


P  =  PCJ  1 
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thus,  when  r  =  r  =  Dt,  p  =  p-,T  according  to  Equation  (47).  The  isentropic 
8  LJ 

expansion  is  terminated  at  the  point  where  the  velocity  u  drops  to  zero 


where  then 


r  *  r2  *  s[aCJ  <y2  +  11  *D(y2‘l)]1 


( 
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p2’pCj[j(y2  +  1)  "  1  (),2  *  l)  i^] 


The  length  of  the  region  occupied  by  the  expansion  is 
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and  increases  in  length  with  time  t.  This  can  also  be  observed  in  Figures 
7,  8,  and  9. 

In  view  of  the  discussion  above  the  variation  of  pressure  at  a  fixed  posi¬ 
tion  along  the  detonation  tube  would  be  as  Indicated  below: 
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p(t)  is  shown  for  two  values  of  r,  r  and  r.  with  r.  >  r  .  For  a  fixed 

a  d  b  a 

radius  the  variation  of  the  pressure  Is  given  by 


(53) 


The  ground  impulse  I  per  unit  area  at  a  fixed  r  is  defined  by  the  integral 


I(r,t)  J  (p  -  Pj)  dt 
JQ 


(54) 


From  Equations  (51),  (52),  and  (53),  it  follows  that 

I(r,  t)  =  0  ;  t<t 
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The  expressions  for  ground  impulse  can  be  reduced  to  the  following  dimen¬ 
sionless  form  in  terms  of  the  similarity  variable  X  =  r/Dt,  and  a  dimension¬ 
less  Impulse  t  given  by 


l(r,t)*0  for  X  >  1 


(58) 
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for 


x2  <  X  <  1 


where 
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For 


X  <  X, 


The  expressions  (58)  through  (60)  are  in  a  universal  form.  The  dimen¬ 
sionless  integrals  only  need  to  be  evaluated  once  for  a  given  y 2  and 
D/a^j.  It  should  be  noted  that  D/a Is  also  equal  to  the  density  ratio 
across  the  C-J  detonation. 

The  dependence  of  I  on  the  radius  of  the  point  where  I  Is  to  be  evalu¬ 
ated  enters  through  the  quantity  t  and  through  the  definition  of  X  as  (r/Dt). 

s 

3.  Chapman  -  Jouguet  Condition 

Evaluation  of  I  from  the  dimensionless  Equations  (58),  (59),  and  (60) 
requires  a  knowledge  of 

D  »  pcA  •  P2^PCJ  *  D^aCJ  *  X2 

and  these  quantities  can  be  determined  by  using  conservation  laws  across 
the  detonation  together  with  the  Jouguet  hypothesis  that  the  velocity  of  the 
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fluid  downstream  of  the  wave  moves  with  sonic  velocity  with  re  sped  to  a 
wave  front.  It  is  here  that  suitable  approximations  can  result  in  great 
simplification. 

For  the  wave  shown  below  (not  necessarily  C-J) 


with  Uy  the  velocity  immediately  downstream  of  the  wave,  the  conservation 
equations  are 

Mass  Pj  D  =  py  (D  -  Uy)  (61) 

2  2 

Momentum  D  +  Pj  *  py  (D  -  Uy)  +  py  (62) 

P,  r>2  PT  (D-uJ2 

Energy  •1  +  ft+T'*y  +  p  + - 2  <63) 

1  y 

where  now  the  velocity  ts  denoted  by  u.  In  Equation  (63),  e^  and  ey  are  the 
specific  Internal  energy  of  the  fluid  including  the  chemical  heats  of  forma¬ 
tion.  If  are  the  reactants  and  are  the  products  in  the  equilibrium 
mixture  downstream  of  the  detonation,  then 
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(64) 
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=  molecular  weight  of  species 

=  number  of  moles  of  each  species  of  reactant  and  products 
=  standard  heat  of  formation 
=  reference  temperature  (e.  g. ,  may  be  0°K) 


Also 
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In  the  case  of  gaseous  detonations,  both  the  reactants  and  products 
satisfy  the  perfect  gas  equation  of  state  although  the  ratio  of  specific  heats, 
the  specific  heats,  and  the  gas  constant  will  change  across  the  wave.  Then 
the  energy  equation  is  often  written  in  the  form 
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and  presumably  Tq  =  0°K  is  used  as  a  reference.  The  quantity  Q  is  related 
to  the  chemical  heat  release  per  unit  mass  of  fuel-oxidizer  mixture  and  is 
given  by 
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It  is  noted  that  Sa^Aj  =  Bt  since  mass  must  be  conserved  in  the 

combustion  reaction. 

The  equations  of  state  on  the  two  sides  of  the  detonation  are  now 
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pi  yi ' 1  1  111  py  yy  - 1  y  y  y  : 


Using  the  perfect  gas  equations  and  the  energy  equation  in  the  form  (63a) 
and  the  Jouguet  Hypothesis  which  requires  that 


s  a  s  a 


CJ 


(70) 


it  is  readily  shown  that  the  density  ratio  across  a  C-J  detonation  is  given  by 


where 

Mp  *  D/a^  ■  Mach  number  of  the  detonation 
yy  *  yCJ  *  *  cofll>t,ult  In  burned  products. 
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o 

For  most  detonations  y^  -  1.  3-1. 4  while  >  3.  Hence,  y^M^  »  1  so 
that 


+  1 


2 


(72) 


i.e. ,  the  density  ratio  is  essentially  independent  of  the  heat  release  across 
the  detonation  or  of  the  detonation  Mach  number  For  many  fuel-air 
mixtures  y^  ~  i.  2  so  that 


PCJ  ->.2.2 

pj  =  1.2 


1.83 


(73) 


In  any  case  Pqj/Pj  does  not  vary  greatly  with  y2  as  is  evident  from  the 
following  tabulation: 


*2 

1.1 

1.2 

1.3 

1.4 

y2  +  l 
y2 

1.91 

1.83 

1.77 

1.71 

Exact  calculation  of  Taylor^  ^  and  Lewis  and  von  Elbe^  yielded  the 
following  results  for  Pgj/py 


Fuel- 

Oxidizer 

H2  +  5°2 

CO  +  O 

C2«2  ♦  10O2 

CH4  +  202 

pcA 

1.78 

1.88 

1.84 

1.904 

For  almost  all  C-J  detonations  it  appears  that  Pqj/Pj  lies  in  the  vicinity  of 
1.8.  In  the  case  of  Hj  +  (1/2)  02,  y2  *  1. 215  so  that  (y2  +  l)/y2  =  1. 82. 
The  approximate  relation  f  Equation  (72)]  thus  appears  to  be  quite  accurate. 

It  is  readily  shown  that  the  pressure  ratio,  P£j/P}>  l0  given  by 


2 

and  Pqj/Pi  will  depend  very  strongly  on  .  Following  the  treatment  in 
Sedov,  it  is  readily  shown  that 

*1  =  ^D2 

satisfies  the  quadratic  equation 


2  A  yl 

Zj  -  2Zjfi  +  — g  =  0 


where 


n2  (v22  - 1) 


6  - (l+a-y. 

Vj  (y,  - 1)  1 

and  where  the  dimensionless  heat  release  per  unit  mass,  ^  is  defined  as 

A _ 9_ 

^  Cn  T1 
P1  1 


The  solution  of  Equation  (75)  is 


and  the  solution  with  the  positive  sign  must  be  rejected  here  since  it  leads 
toMD<l. 

For  almost  all  detonations  Q  »  1.  For  Instance,  assigning  an  approxi¬ 
mate  heating  value  of  19,  000  BTU/lbm  of  fuel  to  methane  yields  Q  -  8  for 
a  stoichiometric  methane  air  mixture  at  520°R.  With  Q  »  1  it  follows 
from  the  definition  of  A  and  Equation  (75)  that 


(71) 


or 


<  v  - » 


%M  2  n  mt 

“d  *  2  (Vj  -  1) 


A 

Q 


D2  s  2(y22  -  1)  Q  (77a) 

A  2 

Again,  since  Q  »  1,  W  »  1,  It  follows  from  Equations  (32)  and  (29) 
that 


(78) 


Actually  when  y..,  y0  and  the  pressure  ratio  are  known,  Equation  (78)  can 

1  1  A 

be  used  to  compute  the  dimensionless  heat  release  Q.  For  Instance,  In 
**»«»  case  of  Hg  +  (1/2)  Og,  =  1.4,  y2  =  1.  215,  P^/Pj  *  18.05  so  that 
s  2.0. 

Finally,  the  temperature  ratio,  T^/T^  across  the  wave  can  be  obtained 
from 


T1  P1  PCJ  R2 


„  (y2_1)  A#t2 

2V2(r1-i)(r2  +  i)Q^ 


(79) 


where  7ft  is  the  molecular  weight. 

4.  Evaluation  of  Ground  Impulse* 

An  evaluation  of  ground  Impulse  will  now  be  considered.  It  Is  necessary 
to  evaluate  the  Integral^ (A): 

1  2y2/(y2-l) 

Jw’(  *  <8o) 

♦Section  C.  4  and  C.  5  are  natural  outgrowths  of  the  foregoing  material. 
The  Information  reported  In  these  two  sections  was  originally  done  for 
Dr.  Norman  Slagg  of  Plcatlnny  Arsenal  under  a  consulting  arrangement. 

For  the  sake  of  continuity  and  completeness,  the  Information  la  also  pre¬ 
sented  here. 


Now  using 


+  1 


2 


J  (X)  becomes 


«#<»•/  Mi+rv<x'i)] 


>2/(y2’1) 


dX 


(80a) 


Although  the  integral  (80a)  can  be  determined  analytically,  it  is  computa¬ 
tionally  simpler  to  evaluate  (80a)  numerically  using,  e.  g. ,  the  trapezoidal 
rule.  The  Integrand  of  (80a)  is  tabulated  for  three  different  values  of  y 
below: 


\ 

j_r1+! 

tr2nr2-V 

2  <> 

A 

A 

v2  '  'J 

y2.1.2 

yj-1.8 

0.5 

1.44 

1.41 

1.390 

0.6 

1.23 

1.22 

1.21 

0.7 

1.11 

1.10 

1.10 

0.8 

1.042 

1.042 

1.040 

0.9 

1.010 

1.010 

1.010 

1.0 

1.000 

1.000 

1.000 

It  can  be  seen  that  over  the  range  of  interest  the  integrand  of  Equation 
(80a)  and  hence  J(X)  is  Independent  of  yg  and  thus  may  be  considered  as  a 
universal  function. 

For  X  <  Xg,  it  is  necessary  to  know  both  Xg  and  p^/p^j.  Using  the 
result  Uqj/D  -  P\/Pqj  *  Vj  /(*  j  ♦  1),  it  follows  that  Xg  -  1/2  so  that  Xg 


is  independent  of  y 2>  Q,  etc.  Also 


CJ 


■[K-K-a  ^**4] 


2V2/(Vj-D 


(81) 


P2/Pcj  ls  tabulated  versus  y2  below 


II 

Cl 

1.1  1.2  1.3 

P2^PCJ  " 

0  .  368  0.  358  0  .  348 

and  it  can  be  seen  that  P2^CJ  *S  es8en^a^y  Independent  of  y 2  in  the  range 
of  interest. 

The  dimension! c  as  impulse  t  now  becomes 
1=0  ,  x> 1 

I  =J(A)  -  (Pl/PcJ)[(l/A)  -  1]  ;  1/2  <  A  <  1  (82) 

I  =^(1/2)  +  0. 358  [(1/X)  -  2J  -  (Pj/Pcj)[(1/A)  -  1)  ;  0  <  X  <  1/2 

Actually  It  is  physically  more  meaningful  to  use  the  similarity  variable 

r  «  1/X  «  t/t 

where  t#  is  tie  time  for  the  wave  to  reach  the  r  where  I  is  to  be  determined. 
Then  a  consideration  of  Equation  (82)  makes  it  clear  that 

U(r)  « t  +  (Pj/PcjMt  -  1)  (83) 

is  independent  of  the  properties  of  the  detonation,  l.  e. ,  of  yv  y2>  C2,  etc. 


U(r)  is  now  considered  to  be  the  universal  impulse  function.  U(r)  is  tabu¬ 
lated  below  and  is  also  plotted  In  Figure  12. 


A 

T 

U(t) 

A 

T 

U(t) 

1.00 

1.00 

0.00 

0.60 

1.667 

0.426 

0.95 

1.050 

0.050 

0.55 

1.818 

0.488 

0.90 

1.110 

0. 1003 

0.50 

2.00 

0.556 

0.85 

1.178 

0.152 

0.40 

2.50 

0.740 

0.80 

1.250 

0.203 

0.30 

3.33 

1.046 

0.75 

1.330 

0.256 

0.20 

5.00 

1.659 

0.70 

1.429 

0.310 

0.10 

10.00 

3.500 

In  terms  of  U(r) 

1  *  "I'.Rf  U(T)  • (T  • l)] 


(84) 


Several  properties  of  I  are  now  readily  determined.  Thus,  since 


s 


r  r  .  PCJ  -  ...  (Yi  '  V 

D  »-«Q  '  pi  Vi-y®  ’ 


and 


Ppj 

U(r)  »  (r  -  1) 


It  follows  that 

I  oc 


■  PCJ  J  Q  Jty^Yt  '  1)(yl  ~ 
pi‘s  P,  "pi r  V  v7  V  ?vT5 


for  a  fixed  r  *  t/tg. 

On  the  other  hand,  It  may  be  desirable  to  know  how  I  Is  related  to  t, 
the  actual  time,  rather  than  the  dimensionless  time  r.  From  Equation 
(82)  It  Is  cl4ar  that  as  t  -  ec  (or  t  »  t^) 


t«(o.  358  -  — V 
V  PCJ' 


Therefore, 

/^CJ\  /  Pi  \  (y2  ‘  ^  a  /  Pi  \ 

1  if  K(° 358  -  ?i)‘  ■  *i  pi(°- 358  -  £)* 

u$»  i, 

(Y2  -  1)  a 

IocO.716^— jyQPjt 


(85) 


(86) 


In  the  case  of  Hg  +  (1/2)  Og 

PCJ/Pl- 18.05  ;  therefore  ,  p^/p^j  =  0.0555 

so  that  neglecting  P^/Pqj  In  Equation  (44)  results  In  an  error  of  the  order 
of  15%  In  the  evaluation  of  I. 

Specific  calculations  of  I  as  a  function  of  time  have  been  made  for 
detonations  In  Hg  +  (1/2)  Og  for  r  «  0,  1,  2,  and  3  meters.  The  results 
are  presented  In  Figure  13  below.  It  can  be  seen  that,  after  an  Initial 
period  representing  the  passage  of  the.  expansion  wave,  the  impulse  varies 
linearly  with  time.  With  Increasing  distance  from  the  center  of  symmetry, 
the  Initial  non-linear  region  broadens  because  of  the  broadening  of  the 
expansion  wave  with  Increasing  time.  The  Initial  slope  of  the  I-t  curves 
Is  steeper  than  the  final  slope  because  Immediately  behind  the  wave  p  =  pCJ 
but  then  drops  to  pg  <  p^j. 

Knowing  p^j,  p^,  and  D  curves  of  I  for  other  detonable  mixtures  are 
readily  computed  using  the  universal  Impulse  function  U(r). 
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20.0 
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5.  Dynamic  Impulse  at  a  Fixed  Point  -  Spherical  Detonations 

As  in  the  planar  case,  a  spherical  detonation  is  followed  by  an  expan¬ 
sion  wave  within  which  the  velocity  of  the  fluid  drops  to  zero  and  the  pres¬ 
sure  drops  from  p^,j  to  some  lower  value  Pg-  The  variation  of  pressure  p 
and  velocity  u  behind  a  spherical  C-J  wave  are  shown  qualitatively  below: 


As  the  wave  moves  outward,  the  width  of  this  expansion  wave  increases. 
Now  the  dynamic  impulse 


I 


d 


(42) 


is  sought  due  to  the  passage  of  a  spherical  C-J  detonation. 

It  can  be  shown  that  for  a  spherical  C-J  detonation  u/u^j,  and  p/p^j 
are  universal  functions  of  r/R  =  r/Dt  =  A  which  depend  only  on  the  value  of 
y g,  the  ratio  of  specific  heats  downstream  of  the  wave.  In  terms  of  the 
dimensionless  variable  A,  the  dynamic  impulse  1^  then  becomes 

1=0  ,  A>1  ,  or  t<t 

d  ’  ’  s 


_  PCJ  °  ts  f  /  p  \  ✓  u  \2  dA 
d  2  J  VPCJ/\UCJ/  A2 

A 

A.  <  A  <  1  or  t  <  t  <  t  (87) 

Cd  S  u 
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At  the  end  of  the  expansion  wave  X  =  Xg  and  then 


0cj  D  *. 


const 


In  the  case  of  spherical  symmetry,  analytical  solutions  for  ( p/p and 
(u/u^j)  are  not  available;  rather,  these  functions  must  be  determined  by 
numerical  integration  of  the  conservation  equations. 

The  variation  of  (u/u^j) ,  (p/Pqj)  >  a/a^j  have  been  computed  by 

Manson  [see  Taylor^]  for  C0H0  +  O,  corresponding  to  y9  =  1. 13  and 
//»\  6  6  6  ® 

by  Sedov'  '  for  y ^  =  1.67.  In  Figures  14  and  15,  p/p^j  and  u/u^j  the 
two  cases  are  compared,  and  it  can  be  seen  that  the  results  are  quite 
insensitive  to  y^.  Hence,  Manson’s  results  will  be  used  to  evaluate  the 
Integrals  in  Equations  (87)  and  (88). 

First,  Equations  (87)  and  (88)  will  be  written  in  a  somewhat  different 
form,  i.e. , 

„  PCJ  .  [  P  /  u  \2  dX 


C-r)  a,  J  J 


PCJ  aCJ  ’  PCJ  (D  '  “CJ*  *  P1  D 


Therefore, 


1 

Vo  +  1 


since 


VPCJ  =  y2/iY2  *  X) 


Hence,  we  can  also  write 


“ 
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since  Ag  =  0.  5. 

The  universal  function  U  (r)  has  been  evaluated  and  Is  shown  in  the 

D 

tabulation  below  and  has  been  plotted  in  Figure  16,  for  spherical  detonations. 


A 

Us 

T 

A 

Us 

r 

1.000 

0 

1.00 

0.900 

0.08836 

1.111 

0.990 

0.01695 

1.010 

0.900 

0.08836 

1.111 

0.980 

0.02965 

1.020 

0.850 

0.1053 

1.178 

0.970 

0.04038 

1.031 

0.80 

0.1183 

1.250 

0.960 

0.04968 

1.041 

0.75 

0.1268 

1.331 

0.950 

0.05795 

1.051 

0.70 

0.1318 

1.430 

0.940 

0.06540 

1.062 

0.65 

0.1345 

1.539 

0.930 

0.07200 

1.073 

0.60 

0.1359 

1.667 

0.920 

0.07793 

1.088 

0.55 

0.1360 

1.819 

0.910 

0.08336 

1.099 

0.50 

0.1362 

2.000 

Id  has  been  computed  for  Hg  +  (1/2)  Og  at  r  =  lm,  2m,  3m  and  is  shown 
in  Figure  17. 

From  Figure  17  it  can  be  seen  that  the  dynamic  Impulse  increases  to 
some  final  value  after  which  it  remains  constant.  This  reflects  the  presence 
of  a  stagnant  region  about  the  center  of  symmetry.  As  r  increases,  the 
width  of  the  expansion  wave  behind  the  detonation  front  increases;  conse¬ 
quently,  the  total  dynamic  Impulse  Increases,  and  this  behavior  is  reflected 
In  the  graph  (Figure  17). 


U  (r)  ■  Universal  Impulse  Function  for  Spherical  Detonation 


r-r/R 

Figure  16.  Universal  Impulse  Function  Us(r)  Spherical  Detonation 


gOi  X3M  to*  'ts|ndui|  D|«miXo  -p| 
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Figure  17.  Variation  of  Dynamic  Impulse  per  Unit  Area 


The  variation  of  the  final  value  of  1^  is  readily  determined  and  is  given 


r(m)  1  2  3  10  20  30  40  100 

^inai  (atm  see  x  W3)  0.  356  0.712  1.068  3.56  7.12  10.68  14.22  35.6 

D.  PROPAGATION  OF  C-J  DETONATIONS  THROUGH  A  CLOUD  OF  DROPS 

In  the  discussion  of  the  C-J  detonation  above,  it  has  been  assumed  that 
both  the  burned  and  unburned  mixture  are  in  gaseous  form.  However,  for 
the  FAE  the  detonation  propagates  into  a  mixture  of  fuel  drop  and  a  gaseous 
oxidizer.  The  combustion  products  downstream  of  the  wave  will  be  in  purely 
gaseous  form.  The  fact  that  the  fuel  drops  do  not  contribute  to  the  pressure 
or  volume  occupied  by  the  unburned  mixture  must  be  considered  in  evaluat¬ 
ing  the  C-J  jump  conditions  across  the  wave.  This  aspect  of  two-phase  det¬ 
onations  has  been  considered  by  Williams^ and  Nicholls  et  a/1^.  The 
Mach  number  of  the  two -phase  detonation  is  lower  and  the  pressure  and 
density  ratios  are  higher  than  for  the  purely  gaseous  detonations. 

In  evaluating  the  spray  detonation  jump  conditions,  it  is  assumed  that 
the  spray  of  density  P|  per  unit  volume  of  mixture  is  uniformly  distributed, 
occupies  negligible  volume,  and  moves  at  the  same  velocity  as  the  gas. 

Then,  If  the  perfect  gas  equation  of  state  is  used  and  the  parameter 


0*1+  (P|/Pi) 


m 


which  Is  the  mass  of  liquid  and  gas  per  unit  mass  of  the  gaseous  component 
is  introduced,  the  equations  of  mass,  momentum,  and  energy  conservation 
across  the  wave  become 

a  (vj  -  c)  *  p2  (v2  -  c)  (93) 

2  2 

pj  +  Of  Pj  (Vj  -  c)  =  p2  +  P2  (v2  -  c)  (94) 


~4(vx-  c)2 


Pi  2 


+  Q 


+ 


(95) 


The  heat  release  Q  is  now  defined  as 

Q  *  $  +  (a  - 1)  e£/a 

with  Q  the  heat  released  per  unit  mass  of  mixture  due  to  chemical  reaction 
and  change  of  phase  and  with  e£  the  specific  internal  energy  of  the  liquid. 

The  self  similar  C-J  solution  described  in  subsection  D,  above,  will 
remain  valid  even  for  two -phase  detonations.  The  presence  of  the  liquid 
fuel  will,  however,  modify  the  values  of  z2  and  V2  and  R2  Immediately 
downstream  of  the  waves. 

E.  THE  INITIATION  PROBLEM 

As  mentioned  above,  the  initiation  of  the  detonation  by  the  secondary 
blast  is  a  crucial  part  of  the  FAE  process.  It  Is  Important  to  know  how 
much  blast  energy  Is  required  to  make  a  fuel  cloud  of  a  given  shape  and 
droplet  density  detonate.  The  answer  to  this  question  will  require  exten¬ 
sive  experimentation;  however,  various  theoretical  aspects  of  the  initiation 
problem  are  discussed  here.  Even  in  the  case  of  gaseous  detonations,  there 
la,  at  present,  no  completely  satisfactory  theory  to  describe  the  blast  initia¬ 
tion  of  a  detonation  wave.  However,  some  of  the  analytical  results  which 
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have  been  obtained  in  the  gaseous  case  should  be  applicable  to  the  FAE 
initiation  problem  and,  therefore,  will  first  be  discussed. 

1.  Blast  Initiation  of  Gaseous  Detonations 

The  problem  of  establishing  what  happens  when  a  blast  wave  is  set  off 
in  a  combustible  gas  mixture  has  received  the  attention  of  a  number  of 
investigators.  The  current  status  of  this  problem  is  discussed  in  the  re¬ 
views  by  Chernyi  et  al^  and  Lee  et  al^,  and  in  the  articles  by  Bishimov 
etal^,  Korobeinikov^,  and  Bach  et  al^ .  When  combustion  reactions 
are  very  rapid,  the  detonation  front  can  be  treated  as  an  exothermic  dis¬ 
continuity.  Then  it  can  be  shown  that  blast  waves  weaken  with  increasing 
radius  and  eventually  change  to  a  C-J  detonation  travelling  at  constant 
speed.  In  the  planar,  v  =  1  case,  this  transition  occurs  only  as  r  —  oc; 
however,  for  cylindrical  and  spherical  geometry  the  transition  occurs  at 
a  finite  value  of  the  radius  of  the  order  of  the  critical  radius  r+.  The 
problem  of  blast  Initiation  when  the  reaction  rates  and  hence  the  detonation 

thickness  are  finite  was  considered  by  Bishimov  et  al^ .  The  results  of 

(2) 

this  analysis,  together  with  experimental  data,  led  Chernyi  et  al  to 
conclude  that: 

M.  .  .  with  the  llmlatlons  of  the  one  dimensional  flow  model 
the  initiation  of  an  overdriven  detonation  wave  by  means  of 
a  piston  (the  role  of  which  may  be  played  by  an  expanding 
compressed  gas) ,  or  by  concentrated  energy  Input,  followed 
by  the  weakening  of  such  waves,  does  not  lead  (in  cases  of 
practical  Interest)  to  their  transition  to  the  Chapman - 
Jouguet  mode. " 

Thus,  with  Increasing  radius  the  width  of  the  reaction  zone  increases  and 
the  combustion  zone  becomes  uncoupled  from  the  leading  shock  wave. 

Regardless  of  these  theoretical  results,  blast  waves  do  initiate  detona¬ 
tions,  and  the  key  to  this  conflict  between  theory  and  experiment  lies  in  the 
portion  of  the  quotation  above  which  has  been  underlined  by  the  present 


author.  Most  detonations  do  not  have  a  one -dimensional  structure; 

rather,  there  are  transverse  waves  within  the  reaction  zone  which  generate 

local  hot  spots  and  thus  continually  renew  the  combustion  process 
(151 

[  Strehlow  ] .  The  blast  Initiation  of  detonations  thus  appears  to  be 

associated  with  instability  of  the  leading  shock  front  and  the  formation  of 
a  complex  three-dimensional  detonation  structure. 


A  phenomenological  theory  which  predicts  the  main  features  of  blast 

(71 

Initiation  has  been  developed  by  Bach  et  al  .  The  key  element  of  this 
theory  is  that  the  interaction  between  the  combustion  and  the  propagation 
of  the  blast  wave  is  accounted  for  by  treating  the  wave  as  an  exothermic 
discontinuity  with  an  effective  heat  release  Q  ,  which  is  only  a  fraction 

Cft 

of  the  heat  of  reaction  Q.  Thus,  Bach  et  al  set 

Qe  =  Q*F(iJ,z)  (96) 


where 


V  = 


» 


(97) 


and  rQ  Is  a  quantity  called  the  explosion  length  which  is  defined  by 


(98) 


ky  *  1,  2t,  4v  for  v  *  1,  2,  3 

In  general,  rQ  is  thus  greater  than  the  critical  radius  r,.  The  choice  of 
the  function  F (fj,  a)  is  empirical  and  depends  on  the  magnitude  of  z  relative 
to  6  *  A/ro,  a  dimensionless  Induction  distance  and  on  the  magnitude  of  V 
relative  to  the  Chapman-Jouguet  value,  and  to  yQ,  the  critical  value 
for  auto^ignition  behind  the  shock.  In  the  theory  the  ignition  delay  or 


induction  distance  A  is  chosen  of  the  order  of  the  ignition  delay  distance 
of  a  plane  C-J  detonation.  To  compute  the  propagation  speed  of  the  wave, 
the  variation  of  velocity  density  and  pressure  behind  the  wave  are  a  priori 
assumed  to  have  the  same  form  as  previous  analytical  blast  wave  solutions. 

A  key  result  of  the  theory  is  that  for  a  large  ignition  energy  Eq  such 
that  6  =  (A/rQ)  ^  1  the  blast  wave  monotonlcally  decays  to  a  C-J  wave. 
When  the  ignition  energy  is  small  so  that  5  »  1,  the  wave  decays  to  an 
acoustic  wave.  Qualitatively,  the  theory  reproduced  the  experimentally 
observed  blast  initiation  of  detonations.  The  critical  energy  Eq  for  initia¬ 
tion  is  thus  related  to  the  ignition  delay  distance  so  that 

(E)  oc  p1a12  P«  £?  (99) 

crit 

(7) 

The  question  now  is  whether  a  theory  similar  to  that  of  Bach  et  al  can 
be  developed  for  the  FAE.  This  problem  is  discussed  qualitatively  below. 

2.  Blast  Initiation  of  Fuel  Droplet  Clouds 

It  is  clear  from  the  discussion  above  that,  regardless  of  the  theoretical 
model  used,  the  ignition  delay  or  induction  distance  A  will  play  a  central 
role.  In  the  gaseous  case  A  is  governed  by  the  chemical  kinetics  of  the 
combustion  reactions;  however,  in  two -phase  combustion  the  process  by 
which  the  liquid  fuel  and  the  gaseous  oxidizer  are  mixed  plays  the  central 
role. 

The  ignition  of  single  fuel  droplets  by  an  Incident  shock  wave  has  been 

fl6)  (17) 

extensively  investigated  by  Kauffman'  ,  Kauffman  et  alv  ,  Kauffman  and 
Nicholls^®\  Ragland^,  Ragland  et  and  Dabora  et  al^^.  After 

the  shock  passes  over  the  droplet  a  detached  bow  shock  is  formed  ahead  of 
the  drop,  the  droplet  flattens  out,  and  a  fine  mlcromtst  is  stripped  from 
the  droplet  by  the  surface  boundary  layer  induced  by  the  flow  over  the  drop. 
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This  micromist  is  entrained  in  the  droplet  wake,  and  the  ignition  process  is 

characterized  by  an  explosion  of  the  fuel  mist  in  the  wake.  Photographs  of 

(171 

fuel  droplets  at  various  stages  of  ignition  are  presented  by  Kauffman  et  al 

Unfortunately,  no  simple  theory  or  correlation  exists  for  the  ignition 
delay  time  or  the  ignition  delay  distance  A.  In  the  case  of  both  reacting 
and  nonreacting  drops  it  is  possible  to  define  a  dimensionless  droplet 
breakup  time 


-  u2  p2 

T  =  —  —  t 

Tb  do  Pl  *b 


(100) 


where  Ug  is  the  velocity  behind  the  shock,  P|  is  the  density  of  the  liquid  fuel, 
dQ  is  the  droplet  diameter,  and  t^  is  the  droplet  breakup  time.  This  dimen¬ 
sionless  breakup  time  was  found  to  be  approximately  constant  over  the 

(171 

range  of  experimental  conditions  investigated  by  Kauffman  et  al  .  For 
dlethycyclohexane  drops  was  found  to  lie  in  the  range  of  3. 0  to  5. 0. 
Further  the  droplet  ignition  delay  time  tras  found  to  be  of  the  same  order 
of  magnitude  as  the  value  of  may  provide  the  basis  for  a  first  estimate 
of  the  ignition  delay  distance  A  for  droplet  clouds. 
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Bach  et  al  '  chose  A  as  the  Ignition  delay  distance  behind  a  steadily 

propagating  C-J  detonation.  If  it  is  supposed  that  t.  ~  t.  ,  then  following 
(71  o  « 

Bach  et  al  As  Dtfe  where  D  is  the  C-J  velocity.  In  the  expression 

[Equation  (100)]  tor  t^,  Ug,  is  the  velocity  induced  by  the  shock  wave 

which  precedes  the  reaction  zone  of  the  detonation  front  as  pictured  below, 

while  Pg  will  be  equal  to  p§,  the  density  behind  the  shock. 


From  Equation  (100)  it  now  follows  that 


A  =  —  T.d  J  — 
u2  bo  V  P2 


From  the  continuity  equation  across  the  shock  wave 


(D/uj)  =[1  -  (P1/P2)] 

so  that  Equation  (101)  can  also  be  written  in  the  form 


A-[i-v<vr1Tb-.V=VS" 


(101) 


(102) 


Combining  Equations  (102)  and  (99)  then  leads  to  the  following  expression 
for  the  critical  initiation  energy: 


(Bo>cri.°CP»irt  1  -  W  J 


(103) 


Thus,  the  critical  Initiation  energy  will,  according  to  Equation  (103),  depend 
upon  Tjy  an  empirical  constant,  upon  the  droplet  diameter  dQ,  upon  the  den¬ 
sity  of  the  gas  ahead  of  the  wave,  on  the  sound  speed  ahead  of  the  wave,  or 
2 

since  a^  =  upon  the  upstream  temperature.  Finally,  (E J 

cr  it 

depends  on  the  shock  density  ratio  Pj/Pj  *nd  on  the  density  of  the  liquid 

fuel.  The  Influence  of  the  detonation  Mach  number,  M  enters  through 

D 


its  influence  on  p^J p^  since  the  shock  density  ratio  is  given  by 


The  validity  of  Equation  (103)  should  be  tested  experimentally. 

There  are  also  certain  geometrical  considerations  which  enter  in  a 
consideration  of  FAE  initiation.  In  the  usual  initiation  analysis  it  is  pre¬ 
sumed  that  the  geometry  of  the  initiating  blast  and  final  detonations  are 
the  same;  however,  this  is  not  necessarily  true.  If  a  point  charge  is  used 
to  detonate  a  cylindrical  fuel  cloud  the  initiating  blast  will  be  spherical, 
as  shown  in  Figure  18  below. 

Transverse  waves  appear  to  play  an  important  role  in  the  initiation  of 

gaseous  detonations  as  discussed  above,  although  a  theory  of  initiation 

which  includes  the  effect  of  these  waves  is  not  yet  available.  The  theory 
(7) 

of  Bach  et  al  only  accounts  for  the  interaction  between  the  combustion 
and  the  initiating  blast  in  a  global  way.  The  detonation  of  two -phase  fuel- 
oxidizer  mixtures  always  appears  to  be  accompanied  by  secondary  shocks 
or  blast  waves  within  the  reaction  zone,  and  an  explosion  or  blast  wave 
appears  to  be  an  essential  feature  of  the  shock  ignition  of  fuel  droplets. 
These  droplet  wake  explosions  should,  in  addition  to  transverse  waves, 
also  play  an  important  role  in  the  initiation  process,  for  the  blast  wave 
emanating  from  the  droplet  wake  provides  a  mechansim  of  interaction 
between  the  droplet  combustion  and  the  initiating  blast.  It  may  be  neces¬ 
sary  to  incorporate  the  effect  of  these  droplet  explosions  in  the  theory  of 
FAE  initiation. 

It  should  be  possible  to  develop  at  least  a  phenomenological  theory  of 
FAE  initiation.  Such  a  theory  will  depend  in  a  crucial  way  on  the  pioper 
description' of  fuel  droplet  ignition. 
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Figure  18.  Initiation  Process  -  Ideal  Versus  Actual 


F.  DEVIATIONS  FROM  THE  IDEAL  CYLINDRICAL  BLAST  MODEL 

Most  of  the  foregoing  analysis  and  discussion  is  based  on  the  idealized 
model  of  the  FAE  consisting  of  a  cylindrical  blast  and  detonation  through  a 
monodisperse  and  uniform  fuel  cloud.  Deviations  from  this  ideal  model 
will  now  be  considered  qualitatively. 


1.  Nonuniform  Distribution  of  Fuel  Drop 

It  is  highly  probable  that  the  blast  dissemination  of  the  liquid  fuel  will 
not  result  in  uniform  cloud  of  fuel  droplets.  Rather,  the  drop  density  and 
hence  the  heat  release  Q  per  unit  mass  of  fuel  oxidizer  mixture  will  be  a 
function  of  radius.  Then  the  self  similar  C-J  wave  discussed  above  will 
no  longer  be  the  limiting  behavior  at  large  time,  i.e. ,  as  t  -  o c. 

However,  self  similar  solutions  are  also  available  when  Q  =  Q(r)  and 
should  be  a  valuable  aid  in  studying  the  effect  of  nonuniform  fuel  distribu¬ 
tion.  For  instance,  suppose  that  Q  varies  with  radius  according  to 


Q  =  B/rw 


(105) 


Then  the  characteristic  parameters  of  the  C-J  wave  problem  are  p^,  p^,  B, 
r,  t,  and  y.  When  0,  there  are  no  combinations  of  the  dimensions  of 
and  Pj  having  the  dimensions  of  the  proportionality  constant  B.  Consequently, 
all  of  the  characteristic  parameters  above  are  dimensionally  independent; 
hence,  no  self  similar  solution  for  the  C-J  detonation  exists,  in  general, 
when  Q  varies  with  r.  However,  If  the  pressure  ratio  across  the  detonation 
is  large  so  that  p^  can  be  neglected,  a  self  similar  solution  can  be  found. 

Then  only  a  single  v  product  can  be  formed  from  the  characteristic  param¬ 
eters  and  leads  to  the  similarity  variable 


X  -  @r/B1/2+"  t2/24" 


(106) 


where  0  is  again  an  arbitrary  constant.  Upon  choosing  the  similarity  vari¬ 
able  X  as  unity  on  the  detonation  front 


r  .  „»/»«>  ,«/*  *•/„ 


(107) 


and 


c  « 


T  ?T3b  /l  p 

.nr an 
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Thus,  r  does  not  Increase  as  rapidly  with  t  than  when  u>  =  0  and  the 
velocity  of  propagation  c  decreases  with  increasing  time. 

A  wide  variety  of  self  similar  solutions  for  variable  Q,  and  p 1  are  also 

(22)  * 
discussed  by  BisMmov  et  alv  ,  and  the  study  of  such  solutions  should  pro¬ 
vide  considerable  insight  to  the  effects  of  nonuniformities  on  FAE  behavior. 

2.  The  Effect  of  Side  Belief 

Even  if  a  uniform  cylindrical  fuel  cloud  can  be  established,  the  upper 
surface  of  the  cloud  will  be  bounded  by  the  inert  atmosphere  rather  than  by 
the  confining  surface  of  the  ideal  model.  As  the  detonation  passes  such  an 
inert  boundary,  a  shock  is  propagated  into  the  inert  gas  and  an  expansion 
usually  propagates  into  the  combustion  products  [Sichel^,  Dabora  et  al^]. 
The  presence  of  the  inert  boundary  causes  the  detonation  wave  to  become 
curved  and  results  in  a  reduction  in  the  velocity  of  propagation. 

A  sectional  view  of  the  possible  form  of  a  cylindrical  wave  with  side 
relief  is  shown  In  Figure  19  below. 


(23^ 

Plane  detonations  with  side  relief  have  been  investigated  by  Sichel  , 

(24)  (25) 

Dabcra  et  al  ,  and  Tsuge  et  al  .  The  author  is  not  aware  of  any 
work  d  "ling  with  the  influence  of  side  relief  on  cylindrical  diverging 
wave  3.  Z1' '  extent  to  which  the  inert  boundary  influences  the  propagation 

of  the  detonation  is  dependent,  to  a  large  degree,  on  the  ratio  of  the  igni¬ 
tion  delay  d< stance  or  reaction  zone  length  A  to  the  width,  h  of  the  explo¬ 
sive,  i.e. .  n  (A/h),  as  well  as  on  the  properties  of  the  inert  boundary. 

An  important  effect  of  the  inert  compressible  boundary  is  that  the  detona¬ 
tion  will  tall  to  propagate  when  (A/h)  exceeds  a  certain  critical  value. 
Clearly,  the  presence  of  inert  compressible  boundaries  will  have  an 
important  effect  on  the  FAE  propagation. 

3.  Propagation  Beyond  the  Fuel  Cloud  Boundary 

Even  in  the  Ideal  cylindrical  model,  a  shock  will  propagate  into  the 
region  beyond  the  fuel  cloud,  1.  e. ,  for  r  >  r  .  This  shock  will  also  con- 
tribute  to  the  total  ground  impulse  provided  by  the  FAE  and  must  be  taken 
Into  account.  As  the  detonation  reaches  the  edge  of  the  fuel  cloud,  a  shock 
will  be  propagated  Into  the  Inert  region  r  >  r  and  either  an  expansion  or 
shock  will  propagate  back  Into  the  burned  gas,  as  shown  in  Figure  20. 
Beyond  the  edge  of  the  fuel  cloud,  r  >  r  ,  the  transmitted  shock  will 

C 

gradually  decay  to  an  acoustic  wave.  It  is  questionable  whether  a  self 

similar  solution  can  be  used  to  describe  the  propagation  of  the  transmitted 

wave  Into  Inert  gas  beyond  the  fuel  cloud.  The  problem  of  the  transmitted 

(26) 

wave  has  been  considered  by  Kiwanv  '  who  used  a  numerical  scheme  to 
compute  the  propagation  of  the  wave  Into  the  region  beyond  the  fuel  cloud. 

4.  Drop  Size  Distribution 

In  the  ideal  FAE  model  it  is  assumed  that  the  droplets  are  uniformly 
distributed  in  the  fuel  cloud  and  that  all  fuel  droplets  have  the  same  diam¬ 
eter.  In  the  actual  fuel  clouds  the  droplet  diameters  will  be  distributed 

, 

over  a  wide  range  of  diameters.  The  smaller  droplets  will  ignl'.e  more 
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Figure  20.  Shock  Interaction  with  Fuel  Cloud  Boundary 

rapidly  than  the  large  droplets  and  no  doubt  will  Influence  the  Ignition  of  the 

larger  droplets.  The  ignition  properties  of  droplets  with  a  blmodal  size 
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distribution  are  currently  being  Investigated  by  Pierce  .  From  his 
results  It  appears  that  the  presence  of  different  sized  droplets  can  have 
an  Important  effect  on  the  Ignition  and  detonation  processes. 

G.  COMPUTATION  OF  C-J  DETONATION  VELOCITIES  FOR  COMPLEX 
HYDRO -CARBON  FUELS 

In  order  to  analyze  FAE  propagation  In  detail,  It  Is  necessary  to  accu¬ 
rately  compute  the  C-J  velocities  for  complex  hydrocarbon -air  mixtures. 
Generally,  such  computations  are  tedious  and  time-consuming.  However, 
the  National  Aeronautics  and  Space  Administration  has,  over  the  past  ten 
years,  developed  a  computer  program  which  very  rapidly  carries  out  such 
computations  for  gaseous  detonations.  This  program  has  been  adapted  to 
the  University  of  Michigan  Computing  Center  and  Is  currently  being  used 
to  compute  detonation  velocities  through  various  MAPP  gas-air  mixtures. 


It  is  planned  to  modify  the  program  to  compute  detonations  through 
droplet  clouds.  The  features  of  the  NASA  program  are  described  in 
detail  below. 

The  NASA  publication  8P-273  describes  the  existing  program 
which  calculates  (1)  complex  chemical  equilibrium  compositions  for 
assigned  thermodynamic  states,  (2)  rocket  performance,  (3)  incident 
and  reflected  shocks,  and  (4)  Chapman -Jouguet  detonations.  It  was  found 
to  be  very  desirable  to  have  parts  (1)  and  (4)  for  use  here. 

The  basic  approach  for  part  (1)  is  as  follows:  Chemical  equilibrium 
is  usually  described  by  either  of  two  equivalent  formulations;  equilibrium 
constants  or  minimization  of  free  energy.  The  program  authors  found  the 
free -energy  minimization  to  be  the  better  of  the  two  in  terms  of  storage 
and  computer  time;  hence  this  approach  is  used. 

The  condition  for  equilibrium  may  be  stated  in  terms  of  any  of  several 
thermodynamic  functions  such  as  the  minimization  of  Gibbs  free  energy 
or  Helmholtz  free  energy  or  the  maximization  of  entropy.  If  one  wishes 
to  use  temperature  and  pressure  to  characterize  the  final  thermodynamic 
state,  the  Gibbs  free  energy  is  most  easily  minimized  inasmuch  as  temper¬ 
ature  and  volume  (or  density). 

The  program  assumes  all  gases  are  ideal  and  that  Interactions  between 
phases  are  neglected. 

The  method  used  for  obtaining  Chapman -Jouguet  detonation  parameters 
uses  three  steps.  The  first  step  consists  of  obtaining  an  initial  estimate  of 
the  detonation  pressure  snd  temperature.  The  second  step  is  an  improved 
estimate  of  these  parameters  by  means  of  a  recursion  formula.  The  third 
step  consists  of  obtaining  the  correct  values  by  means  of  a  Newton -Raphson 
Iteration  procedure.  The  standard  conservation  equations  for  continuity 


£ 

momentum  and  energy  apply  with  the  additional  constraint  that  the  burned 

gas  velocity  Is  sonic  In  the  Chapman -Jouguet  plane. 

1 

The  program  had  to  be  modified  since  It  was  originally  written  for  use 
on  an  IBM  7090  computer  under  a  standard  IBM  operating  system.  [  The 
University  of  Michigan  uses  an  operating  system  called  MTS  (Michigan 
Terminal  System)  and  IBM  360/67  computers.  ]  The  modifications  to  the 
program  have  been  completed,  and  the  C-J  velocities  In  MAPP-air  mixtures 
are  currently  being  computed. 

H.  CONCLUSIONS— THEORETICAL  ANALYSIS 

Initially  the  FAE  behaves  like  a  strong  blast  wave;  then  after  a  period 
of  transition,  the  FAE  tends  toward  the  behavior  of  a  C-J  wave.  Although 
a  completely  detailed  description  of  the  FAE  will  require  lengthy  numerical 
integration  of  the  conservation  equations,  the  discussion  above  indicates 
that  much  useful  Information  can  be  gained,  with  minimal  computational 
effort,  from  a  consideration  of  V  e  well-known  self  similar  blast  wave  and 
C-J  solutions.  For  many  considerations  it  may  be  sufficient  to  model  fuel 
air  explosions  by  a  self  similar  blast  wave  for  r  <  r*  followed  by  a  C-J 
wave  for  rg  >  r#. 

The  ignition  delay  distance  plays  a  key  role  both  in  the  initiation 
process  and  in  establishing  the  influence  of  side  relief.  In  the  FAE  this 
distance  will  be  mainly  governed  by  the  fuel  droplet  breakup  and  ignition 
process  rather  than  by  reaction  rates  as  in  the  gaseous  case. 
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SECTION  m 

EXPERIMENTAL  RESEARCH 

A.  INTRODUCTION 

The  purpose  of  this  phase  of  research  has  been  to  conduct  controlled 
experiments  on  the  cylindrical  propagation  of  two -phase  detonation. 
Heretofore,  controlled  experiments  have  never  been  conducted  on  the 
cylindrical  or  spherical  propagation  of  such  explosions  wherein  fuel -air 
mixture  ratio  and/or  drop  sizes  are  controlled.  Since  such  conditions 
are  Important  in  fuel -air  explosions,  an  estimate  of  their  Influence  is 
desirable.  In  order  to  increase  the  confidence  of  the  theoretical  estimates, 
experiments  were  performed  on  a  scaled-down  model  of  a  fuel -air  explo¬ 
sion  in  the  laboratory.  A  cylindrically  shaped  cloud  was  modeled  which 
is  idealized  to  the  extent  that  cloud  conditions  (l.  e. ,  mean  drop  size  and 
local  mixture  ratio)  vary  only  with  radius  rather  than  with  azimuth  or  in  a 
direction  parallel  to  the  axis.  Due  to  the  axial  symmetry  of  the  problem 
and  in  order  to  minimize  the  magnitude  of  the  explosion,  the  test  chamber 
configuration  chosen  is  a  pie -shaped  segment  of  a  cylindrical  cloud 
(Figure  21). 

B.  DESIGN  OF  THE  EXPERIMENTAL  FACILITY 

1.  Selection  of  Teat  Chamber  Dimensions  and  Fuel  Drop  Size 

The  walls  of  a  detonation  test  chamber  will,  in  general,  have  some 
effect  on  the  detonation  process  due  to  heat  and  momentum  losses  to  the 
wall.  While  a  very  large  chamber  will  essentially  eliminate  the  wall 
effects  on  the  overall  process,  if  a  test  chamber  is  to  be  operated  within 
a  laboratory,  It  Is  Important  that  the  chamber  not  be  excessively  large. 

(21) 

Two-phase  detonation  studies  in  constant  area  test  chambers'  '  have 
demonstrated  that  there  is  a  relation  between  fuel  drop  size,  test  chamber 
dimensions,  and  the  magnitude  of  the  wall  effects  on  the  detonation  process. 


Therefore,  in  the  process  of  determining  the  dimensions  of  a  test  chamber 
which  will  minimize  wall  effects  and  still  be  of  a  practical  size  and  shape, 
it  is  necessary  to  simultaneously  consider  the  size  of  the  fuel  drops. 

The  selection  of  the  size  and  distribution  of  the  fuel  drops  within  the 
test  chamber  must  also  be  based  on  the  practical  limitations  of  the  drop- 
producing  mechanism  and  the  fuel -to -air  ratio  desired. 

The  mass  of  liquid  fuel  per  unit  volume  of  air,  m^,  is  given  by 

mjL  =  N(»/6)  D3  p£  (109) 

where  N  =  number  of  drops  per  unit  volume  of  air 
D  =  drop  diameter 
Pj  *  density  of  liquid  fuel 

The  mass  of  the  air  per  unit  volume  may  then  be  given  by 

ra#  *  [  1  -  N<»/6)  D3]  p#  (110) 

Neglecting  the  volume  of  the  droplet  compared  to  that  of  the  control  volume 
and  neglecting  evaporation  of  fuel,  the  fuel/air  mixture  ratio  from  Equations 
(109)  and  (110)  becomes 

m,/ma  -  N(»/«)  D3  pt/pt  (111) 

The  mixture  ratio,  by  weight,  is  equivalent  to  the  product  of  the  fuel/air 
mole  fraction,  q,  times  the  ratio  of  their  molecular  weights, 

This  Is  expressed  as 

"a,/®,  -  -  n+(l,/Aa> 

Here  n  Is  the  stoichiometric  mole  fraction  and  +  the  equivalence  ratio. 
Equating  (111)  and  (119  and  solving  for  N  gives 

•  (u3> 


The  fuel  selected  for  this  experimental  work  was  kerosene  because  of 
its  relatively  low  vapor  pressure,  reasonably  good  other  physical  character¬ 
istics,  and  its  availability.  For  kerosene  drops  in  theoretical  air  at 
standard  conditions,  we  have  from  Equation  (113)  the  result 

N  =  1. 98  x  10”4  (^/D3)  (114) 

where  drop  size  D  is  given  in  centimeters. 

The  distribution  of  fuel  drops  could  be  rather  simply  determined  by 
Equation  (114)  if  the  drops  were  equally  spread  in  all  directions.  However, 
with  the  hypodermic  needle  technique  of  producing  fuel  drops,  which  was 
contemplated  for  this  work,  the  vertical  distance  between  the  drops  is 
necessarily  relatively  short.  It  is  necessary  for  the  horizontal  spacing 
to  be  appreciably  greater  than  the  vertical  spacing  if  very  rich  fuel/air 
mixtures  are  to  be  avoided. 


The  rectangular  volume  element  which  surrounds  and  is  apportioned 
to  each  fuel  drop  is  defined  as  the  product  of  the  sides  of  the  element. 

Thus 

=  V  *  volume  of  element/fuel  drop  (115) 

x  y  z 

and  V  =  1/N  . 

The  arrangement  of  a  volume  element  within  the  test  chamber  configuration 

planned  is  shown  by  Figure  22.  The  magnitude  of  S  ,  S  ,  and  S  are  not, 

x  y  z 

in  general,  equal,  but  the  following  relations  can  be  assumed: 


(116) 


Using  Equation  (115),  the  number  of  drope/unlt  volume  can  be  expressed  by 


» 


(111) 


i 


In  order  to  arrive  at  practical  values  for  a  and  b  in  Equation  (117),  the 
mechanism  of  drop  formation  must  be  considered.  The  drop -generating 
technique  under  consideration  involves  the  breakup  of  a  stream  of  liquid 
fuel  from  a  small  tube.  Although  such  small  jets  will  generally  break  up 
due  to  instabilities,  the  breakup  would  normally  be  rather  Irregular  unless 
special  conditions  exist.  If  the  liquid  jet  is  disturbed  at  the  frequency  at 
which  it  is  most  unstable,  the  Rayleigh  disturbance  frequency,  the  jet  will 
break  up  and  form  quite  uniformly  sized  drops.  At  such  a  disturbance 
frequency  the  wavelength,  X,  is  given  as  a  function  of  drop  diameter  by 

X/D  =  2.  38  (118) 

Tests  show  X  compares  approximately  with  the  vertical  drop  spacing,  S 

z 

if,  as  assumed  here,  the  terminal  drop  velocity  is  equal  to  the  liquid  jet 
velocity.  The  drop  size  may  then  be  related  to  vertical  drop  spacing 
through 

S  /D  »  2.  38  (119) 

z 

Using  Equation  (117)  in  Equation  (119),  the  number  of  drops  per  unit  volume 
is  related  to  drop  size  by 

N  *»  1/ab  (2.  38fi)3  (120) 

Equating  the  results  of  the  calculations  for  N,  Equation  (114),  found  by 
relating  mixture  ratio  expressions,  and  Equation  (120),  found  on  the  basis 
of  optimum  drop  formation,  the  result  is 

ab  =  375. 0/4  .  (121) 

For  a  stoichiometric  reaction  (4  *  1)  the  product  ab  Is  given  as 

ab*  375.0  .  (122) 


It  follows  that  the  drop  spacing  dimensions  of  Equation  (116)  may  be 
arrived  at  by  selecting  a  desired  drop  size  and  choosing  values  for  the 
constants  a  and  b,  subject  to  Equation  (122).  For  a  given  drop  diameter, 
the  vertical  dimension  is  found  using  Equation  (119);  then  using  the  result, 
Equations  (122)  and  (116),  the  x  and  y  dimensions  can  easily  be  found. 

The  nominal  drop  size  was  selected  following  an  examination  of  the 
frictional  and  heat  losses  to  the  walls  of  a  prospective  test  chamber  iitd 
the  reaction  zone  length  of  a  detonation  process.  The  wall  losses  may  be 
taken  Into  account  through  the  use  of  the  Equation  (123)^^,  where  ug  is 
the  actual  detonation  velocity,  ugQ  the  ideal  C-J  velocity,  y^  the  ratio  of 
specific  heats  at  the  Chapman-Jouguet  plane,  pj  the  fuel  density,  the 
oxidizer  density  In  front  of  the  shock  wave,  the  heat  transfer  coefficient, 
and  r^  the  hydraulic  radius. 

2  1/2  '1/2 

"g/Ugo -[1  +  (SOD/r^  y*  fy/p/'2  C„]  (123) 


The  product  in  the  brackets  Is  small  compared  to  one,  allowing  (u  -  u  1/ 

SO  D 


u_v  the  dimensionless  velocity  deficit  to  be  expressed  by  the  binomial 

1/2 


so’ 
expansion 


Au  /u 
s'  so 


■  (us/u  J  *  15D/rh  O  CH  WpJ 


<12« 


•3 

For  yg  =  1.2,  CH  =  2. 5  x  10  and  the  appropriate  densities,  Equation  (124) 
becomes 


AUs/Uso  *  5-  5  x  10’5  <D/rh>  (125) 

for  D  In  microns  and  rfa  In  inches.  Equation  (124)  is  plotted  in  Figure  23. 
Hydraulic  radius,  defined  as  2x  crossectlonal  area/perimeter,  may  be 
expressed  as  a  function  of  test  chamber  geometry  through  (see  Figure  24) , 

rh  *  2  tan  (6/2) /[(l/r)  +  (2/h)  tan  (6/2)]  (126) 
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Percent  Velocity  Deficit 


The  velocity  deficit  of  Equation  (125)  using  Equation  (126)  can  then  be 
expressed  as  a  function  o i  drop  diameter  and  test  chamber  geometry.  The 
resulting  expression  was  examined,  and  plots  similar  to  Figure  25  were 
made  for  various  chamber  configurations.  For  a  selected  chamber  geom¬ 
etry,  this  figure  gives  chamber  radial  distance  versus  drop  diameter  for 
curves  of  constant  velocity  deficit.  As  expected,  the  velocity  deficit  be¬ 
comes  very  large  in  the  chamber  vertex  region.  It  was  decided  that 
conditions  should  be  selected  which  would  allow  the  flow  properties  to  be 
well  established,  with  few  wall  losses,  by  a  radial  distance  of  6  inches. 
Hence,  a  practical  chamber  geometry  was  selected  on  the  basis  of  a 
tolerable  velocity  deficit. 


Figure  24.  Reference  need  to  Expreee  Hydraulic  Radius  in 
Terms  of  Chamber  Geometry 


4* 


The  nominal  drop  size  was  established  by  considering  the  detonation 

reaction  zone  lengths  as  a  function  of  drop  size  in  conjunction  with  the 

20 

above  mentioned  figures.  It  had  previously  been  concluded  that  the 
reaction  zone  length  for  two -phase  detonations  is  controlled  by  the  drop 
breakup  time.  Breakup  time,  t^,  of  inert  drops  by  shock  waves  can  be 
related  to  the  dynamic  pressure  of  the  convective  flow  behind  the  shock  and 
the  Initial  drop  diameter  D  as  follows^ 

1/2 

t/D-kfp/tpjUj2)]  (127) 

where  k  is  approximately  equal  to  5  and  Pj  is  the  liquid  density,  and  P2  and 
Ug  are  the  convective  flow  density  and  velocity.  A  better  estimate  of 
actual  drop  breakup  time  In  a  detonation  is  achieved  if  an  average  dynamic 
pressure  is  used.  Since  dynamic  pressure  in  a  convective  flow  behind  a  C-J 
plane  is  less  than  10  per  cent  of  that  behind  a  shock  traveling  at  the  same 
Mach  number,  for  M  >  3  an  average  dynamic  pressure  equal  to  one -half 
that  used  in  Equation  (127),  can  then  be  used.  Then  Equation  (127)  becomes 

t/D  -  1. 41k  (p/p,) 1/2  (p,/pj)1/2  (l/Uj)  .  (128) 

In  terms  of  shock  velocity  we  have 

*b*  *b“a/D  *  1,41k  <Vpl)t/2  (pl/p2>1/2  (l  ’  *  (129) 

Equation  (129)  gives  2^  *  162  for  M  «=  3  and  decreases  to  120  for  M  becoming 
Infinitely  large.  For  this  case  then  Equation  (129)  becomes 

3^  ■  tbuf/D  =  a/D  a  162  .  (130) 

In  order  to  have  a  smoothly  propagating  detonation  wave,  the  drop  spacing 
8^  (Figure  22)  must  be  of  the  order  of  the  reaction  zone  length  x  given  In 
Equation  (130).  Thus,  Equation  (130)  becomes 


Sy  <  162D 


(131) 


The  final  selection  of  the  nominal  fuel  drop  size  and  the  test  chamber  dimen¬ 
sions  were  made  on  the  basis  of  allowable  limits  on  velocity  deficit  and  the 
practical  problems  of  drop  generation.  By  examining  a  series  of  graphs 
similar  to  Figure  25  and  by  considering  that  the  allowable  limit  on  the  ve¬ 
locity  deficit  at  a  six -inch  radial  distance  should  be  about  5  per  cent,  a 
nominal  drop  size  of  less  than  400  microns  was  obtained.  By  orienting  the 
pie-segment  shaped  test  chamber  as  shown  in  Figure  22  and  injecting  the 
fuel  drops  through  the  top  surface,  the  vertical  separation  between  the 

drops  Is  S  .  As  a  practical  matter,  the  needle  positions  were  chosen  so 
z 

that  groups  of  fuel  injection  needles  would  be  installed  in  the  top  of  the 
chamber  along  lines  parallel  to  the  x  axis  (Figure  22).  The  separation 
between  the  individual  needles  in  each  of  these  rows  is  Sx  while  the  (pro¬ 
jected)  separation  between  each  row  of  needles  is  S^. 

The  value  of  Sw  is  essentially  fixed  by  the  choice  of  drop  size,  accord- 
z 

ing  to  Equation  (119).  The  value  of  Sy  is  limited  by  Equation  (131).  Com¬ 
bining  values  of  S„  and  S  with  Equations  (116)  and  (122)  results  in  a  value 
z  y 

for  S^.  The  drop  spaclngs  finally  arrived  at  were 

«  0. 64  cm  ;  S  =  5.  4  cm  ;  S  =  0. 095  cm 
x  y  z 

Since  the  chamber  is  essentially  2  Inches  wide,  an  x-separatlon  of  0. 64  cm 
results  in  seven  needles  per  row.  It  was  decided  that  an  overall  chamber 
length  of  about  SO  inches  would  be  suitable  for  this  research.  Therefore, 
eleven  rows  of  needles  with  a  y -separation  of  5. 4  cm  were  required,  or  a 
total  of  77  needles  is  required.  It  is  of  interest  to  note  that  for  a  nominal 
drop  slss  of  400  microns,  the  approximate  length  of  the  reaction  zone  is 
6. 35  cm.  Figure  20  Is  provided  which  shows  velocity  deficit  as  a  function 
of  radial  distance  for  the  selected  chamber  geometry. 


Hydraulic  Radius 


2.  Method  of  Fuel  Cloud  Production 


This  phase  of  the  research  has  been  directed  toward  establishing  a 

satisfactory  technique  for  generating  a  cloud  of  fuel  drops,  initially  of 

uniform  size.  A  drop -generation  technique,  different  from  those  previously 

employed  in  our  work,  was  developed  to  produce  a  spray  of  fuel  drops  of 

nearly  uniform  size  and  spacing  (Figure  29).  The  technique  provides  for  a 

further  extension  to  the  vibrating  capillary  jet  methods,  based  essentially 
(29) 

upon  Rayleigh's  analysis  of  jet  instability.  The  decision  to  develop  a 
new  technique  was  made  following  preliminary  tests  involving  those  already 
in  existence. 

The  basis  of  all  the  techniques  is  to  apply  a  sinusoidal  disturbance  to  a 
liquid  jet  emerging  from  a  vertically  mounted  capillary  tube.  In  each  such 
jet,  the  disturbance  will  grow  in  the  fluid  as  it  moves  away  from  the  capil¬ 
lary  exit.  If  this  growth  rate  is  sufficient,  the  jet  will  break  up  at  its  con¬ 
tractions,  forming  separate  drops.  Rayleigh's  theory  predicts  maximum 
amplification  of  cyclic  disturbances  will  occur  at  those  frequencies  (referred 
to  as  Rayleigh  frequencies,  fR)  determined  by  the  relation  Involving  jet 
velocity,  Uj  and  jet  diameter  d^: 

=  Uj/4.  508  dj  (132) 

Formation  of  drops  of  uniform  size  and  spacing  occurs  more  readily  at  such 
Rayleigh  frequencies.  Furthermore,  at  this  optimum  frequency  the  capil¬ 
lary  tube  inside  diameter  and  the  resulting  drop  size  are  related  through 

D/dj  *  1. 89  (133 ) 

Similarly,  the  wavelength  associated  with  f^  In  Equation  (133)  is  related  to 
drop  size  by 

(134) 


A/D  «  2. 38 


It  should  be  pointed  out  that  Rayleigh's  analysis  was  developed  for  a  station¬ 
ary  capillary  column  of  an  inviscid  fluid.  The  analysis  was  later  extended 
to  include  the  effects  of  liquid  viscosity  and  jet  velocity.  Two  character¬ 
istic  velocities  can  be  associated  with  the  problem  of  liquid  exiting  from  a 
capillary  tube.  The  velocity  of  a  smooth  jet  exiting  from  a  capillary  tube 
is  dependent  on  surface  tension  a,  density  of  the  liquid  P|,  and  diameter 
of  the  jet  dj  through 

Uj2  >  8a/pdj  (135) 

Hence,  Equation  (135)  expresses  the  fact  that  a  capillary  jet  must  have  a 
certain  minimum  velocity,  u  to  avoid  dripping.  Using  Equations  (132) 
and  (135),  a  Rayleigh  frequency  associated  with  this  minimum  velocity 
can  be  obtained 


^-0.627  [<r/(pd^] 


(136) 


Upon  exiting  from  the  capillary  tube,  the  liquid  jet  will  break  up  and  the 
drops  will  accelerate  due  to  gravity  and  approach  their  terminal  velocity. 

In  the  approximate  range  of  drop  diameters  290 p  <  D  <  2600 ^  the  drag 

law  for  falling  drops  (10  <  Re  <  1000)  can  be  taken  as  Re1^2  C-.  *  9. 64. 

(30)  u 

Thus,  it  can  be  shown  that  when  the  initial  jet  velocity  equals  the  termi¬ 
nal  drop  velocity,  the  particular  Rayleigh  frequency  for  these  conditions, 
fp*,  is  determined  by 

ij*  .  0. 112  (p/gV(pg^)l/S  .  (137) 

Operating  at  fp*  tends  to  deter  coalescence  of  successive  kerosene  drops 
greater  than  approximately  350  microns,  since  the  drops  are  formed  at 
their  terminal  velocity  and  therefore  do  not  accelerate  further.  Figure  27 
compares  the  terminal  velocities  of  liquid  drops  In  air,  Uj,  with  the  mini¬ 
mum  jet  velocities  required  to  avoid  dripping. 


DROP  DIAMETER  (microns) 


Figure  17.  Terminal  and  Minimum  Velocities  of  Water, 
Kerosene  and  DECH  Drops 


Heretofore,  the  research  efforts  of  these  laboratories  into  two-phase 

detonations  have  made  extensive  use  of  a  method  of  spray  production 

(31  321 

developed  by  Dabora  et  al  ’  ' .  In  his  method  the  cyclical  disturbance 

is  introduced  in  the  liquid  supply  reservoir  by  means  of  a  vibrating  dia¬ 
phragm.  This  disturbance  is  transmitted  through  the  fluid  moving  in  the 
capillary  needles,  located  directly  below  the  diaphragm,  to  the  free  jet. 

In  his  research  Dabora  chose  to  operate  his  system  at  the  minimum  jei 
velocity.  Consequently,  he  had  drop  coalescence  occurring  for  drop 
sizes  with  umln  <  u.  (Figure  27).  This  was  partially  overcome  in 

such  cases  by  Introducing  an  air  stream  collnear  to  the  liquid  jet,  termed 

(21 

coflow.  More  recently,  T.  H.  Pierce  ’  developed  a  second  method  during 
the  course  of  his  work  involving  both  monodisperse  and  polydlsperse 
spray s^^.  By  his  technique  the  fluid  flow  to  each  individual  needle  is 
carried  by  a  small  inside  diameter  flexible  plastic  tube  which  passes 
through  a  vibrating  platform.  The  oscillatory  disturbance  Introduced  into 
the  fluid  In  this  way  is  transmitted  to  the  jet  Issuing  from  the  needle. 

The  operating  frequency  of  the  disturbance  Is  established  by  Equation  (137), 
l.  e. ,  the  maximum  amplification  frequency  is  chosen  as  the  operating 
frequency,  f  *  f^*. 

It  was  anticipated  from  the  outset  that  the  present  research  would  be 
Involved  with  capillary  needles  numbering  nearly  an  order  of  magnitude 
more  than  had  previously  been  used.  Hence,  every  possible  simplifica¬ 
tion  was  deemed  desirable  as  prospective  spray  generation  techniques 
were  considered.  The  latter  method  by  Pierce  was  rejected  in  view  of 
the  complexity  of  the  plumbing  required.  The  former  method  by  Dabora 
et  al  had  the  undesirable  trait  of  requiring  coflow  to  establish  drops  of 
uniform  size  and  spacing  in  the  present  research.  Undue  eddies  and 
turbulence  would  have  been  established  within  the  test  chamber,  the 


effects  of  which  would  have  been  difficult  to  ascertain.  A  further  word  on 
this  subject  is  in  order.  In  previous  applications  the  direction  of  the  coflow 
and  the  direction  of  detonation  wave  propagation  were  colinear.  In  the 
present  research,  however,  had  coflow  been  used  it  would  have  been  nor¬ 
mal  to  the  wave  propagation  direction.  Hence,  coflow -induced  turbulence 
greatly  enhanced  spray  distribution  uniformity  in  previous  research,  where¬ 
as  in  the  present  work  it  would  have  greatly  disturbed  the  distribution. 

Both  of  these  drop -generation  techniques  make  use  of  a  vibrator  to 
supply  the  necessary  sinusoidal  disturbance.  In  each  case  the  mass  of  fluid 
being  disturbed  is  small;  hence,  the  vibrator  power  requirement  was  never 
questioned  since  a  suitable  vibrator  was  always  available.  However,  in  the 
present  research,  the  mass  of  the  fluid  to  be  disturbed  is  of  an  order  of 
magnitude  greater,  suggesting  larger  vibrator  power  requirements.  Indeed, 
tests  proved  this  statement  to  be  true;  sufficient  disturbances  could  not  be 
obtained  using  in-house  equipment.  Hence,  in  addition  to  the  undesirable 
traits  mentioned,  using  either  of  these  drop -generation  techniques  in  the 
present  research  would  have  carried  with  it  a  demand  for  very  high  vibrator 
power. 

In  view  of  the  foregoing  considerations,  an  alternate  technique  was 
devised  to  generate  the  required  monodisperse  spray.  A  motor -driven 
system  was  developed  in  order  to  generate  stronger  pressure  pulses. 

The  sinusoidal  disturbances  of  the  liquid  fuel  is  generated  by  periodically 
interrupting  the  flow  to  the  needles.  By  this  new  technique,  the  disturbance 
is  introduced  to  the  fluid  through  rotation  of  a  48 -tooth  gear  enclosed  In  a 
cylindrical  case.  The  Inlet  and  outlet  ports  are  so  positioned  that  the  top 
land  of  the  rotating  gear  simultaneously  opens  or  closes  both  ports.  Kero¬ 
sene  under  pressure  Is  fed  through  one  port  to  the  rotating  gear.  The  fuel 
leaving  the  pulse  generator  enters  the  test  chamber  through  a  connected 
series  of  eleven  sub-manifolds,  made  up  of  seven  capillary  needles  each. 


Figures  31  through  34  are  descriptive  of  the  new  drop -generation  system. 
The  new  system  underwent  thorough  testing  to  determine  operating  condi¬ 
tions  necessary  to  generate  drops  of  uniform  size  and  spacing.  Table  I 
compares  the  values  of  drop  size,  operating  frequency,  and  fuel/air  mix¬ 
ture  ratio  obtained  experimentally  with  those  computed  by  Equations  (133) 
and  (137). 

TABLE  I.  COMPARISON  OF  EXPERIMENTAL  DROP  FORMATION 
RESULTS  WITH  THEORETICAL  PREDICTIONS 

Per  cent 


D 

n 


Experimental 
360 fx 
1584  cps 
1/15.4 


Theoretical 
384m 
1550  cps 
1/15.0 

(stoichiometric) 


Difference 

-6.25 

+1.94 

-2.6 


Strain  gage  pressure  transducers  were  located  in  the  system  to  measure 
both  the  static  pressure  of  the  fuel  and  the  magnitude  of  the  pressure  dis¬ 
turbances.  The  pressure  fluctuations  at  the  capillary  needles  were  found 
to  have  a  10  per  cent  peak -peak  variation  about  a  static  pressure  of 
13.  5  pslg. 

Photographs  were  taken  of  the  drops  formed  at  each  of  the  eleven  sub¬ 
manifold  stations  along  the  main  fuel  manifold  to  establish  the  presence  of 
any  nodal  points  in  the  main  manifold  exhibited  by  the  formation  of  drops 
of  highly  singular  sizes  and  spacings.  Such  tests  revealed  the  presence 
of  no  such  points.  Drop  formations  at  the  given  operating  conditions  at 
all  eleven  stations  were  very  satisfactory.  Figure  28  shows  the  type  of 
drop  formation  resulting  when  no  disturbances  are  present,  while  Figure 
29  reveals  the  nature  of  good  drop  formations. 

Inherent  advantages  of  this  drop -generation  system  lie  in  its  ease  of 
maintenance  and  flexibility  In  terms  of  future  test  possibilities.  The 
plugging  of  an  all  capillary  needles  by  combustion  products  has  always 


Figure  28.  Drops  Generated  by  One  Sub-manifold  with  No 
Pressure  Disturbances  Present 
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Figure  29.  Drops  Generated  by  One  Sub-manifold  (Figure  34) 
of  a  Completely  Integrated  System  with  Pressure  Disturbances 
Present.  (Kerosene,  f  =  1580  cps,  360  micron  drops,  1/16  OD 

Capillary  Needles). 
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been  a  problem.  The  use  of  quick  release  fittings  to  hold  the  sub-manifolds 
in  place  allows  for  their  easy  replacement.  With  such  a  capability,  it  is 
also  possible  to  interchange  currently  used  sub-manifolds  with  others  of 
varying  design.  Such  design  changes  could  include  reduction  in  the  number 
of  needles  per  sub-manifold  or  the  use  of  capillary  tubing  of  varying  inside 
diameter,  with  outside  diameter  held  constant.  The  second  change  would 
result  in  tests  involving  varying  drop  size,  while  both  would  involve  tests 
investigating  the  effects  of  varying  mixture  ratios. 

3.  Explosion  Initiation 

The  initiation  of  detonation  in  gaseous  fuel-oxygen  mixtures  is  relatively 
easy  and  many  methods  have  been  used  (l.  e. ,  spark  plugs,  glow  plugs,  shock 
waves,  and  even  cigars).  The  replacement  of  oxygen  with  air  renders  initia¬ 
tion  more  difficult,  yet  many  detonation  experiments  on  gaseous  fuel -air 
mixtures  have  been  conducted.  When  the  fuel  is  in  the  form  of  liquid  drops 
rather  than  gaseous,  initiation  is  apparently  more  difficult.  In  recent  studies 
on  two -phase  detonation,  wherein  the  oxidizer  was  systematically  varied 
from  pure  oxygen  to  oxygen -nitrogen  mixtures,  it  was  found  extremely  diffi¬ 
cult  to  ignite  the  fuel  when  the  volumetric  percentage  of  oxygen  in  the  oxi¬ 
dizer  was  less  than  about  35  per  cent.  Further,  IIT  investigators  were 
unable  to  generate  a  sustained  detonation  with  liquid  fuel  drops  of  several 
types  in  air,  although  it  is  possible  that  part  of  this  difficulty  could  be 
attributed  to  small  cloud  size  and  insufficient  igniter  energy.  In  view  of 
the  above,  it  was  deemed  Important  to  carefully  consider  the  method  of 
initiation  in  connection  with  the  present  research. 

Initiation  techniques  considered  Included  auxiliary  shock  or  detonation 
tubes  as  well  as  solid  explosives.  Auxiliary  shock  and  detonation  tubes 
require  more  support  equipment  than  the  comparable  solid  explosive  sys¬ 
tem,  and  the  time  between  experimental  runs  would  be  longer.  An  advan¬ 
tage,  however,  would  be  that  these  techniques  can  produce  a  sustained 


high  pressure,  high  temperature,  and  high  velocity  blast  directed  at  the 
drops.  On  the  other  hand,  the  explosive  process  is  accompanied  by  a 
rapid  decrease  in  these  gas  dynamic  characteristics  immediately  behind 
the  shock.  The  latter  would  render  initiation  more  difficult.  In  spite  of 
the  latter,  it  was  decided  to  select  the  solid  explosive  technique  due  to 
its  simplicity,  flexibility,  and  the  ease  with  which  the  energy  level  could 
be  increased.  Uncertainty  in  the  location  of  the  energy  threshold  point 
required  experimental  runs  of  vastly  differing  energy  releases.  Blasting 
caps  and  Dupont  Detasheet  C  were  used  in  varying  combinations  to  achieve 
varying  amounts  of  total  energy  release.  This  energy  release  can  be 
established.  The  composition  of  Detasheet  C  is  63  per  cent  Pentaerythrite 
Tetranitrate  (PETN),  8  per  cent  Nitrocellulose  (NQ>  and  25  per  cent  non* 
combustibles.  The  heats  of  combustion  '  for  PETN  and  NC  are,  respec¬ 
tively,  4280  ft -lb/gram  and  2640  ft-lb/gram.  Taking  a  weighted  average 
of  these,  the  heat  of  combustion  for  Detasheet  is  given  as  2911  ft-lb/gram 
or  188  ft  -lb/grain.  Detasheet  was  exploded  using  Dupont  blasting  caps 
with  a  nominal  energy  release  of  1110  ft -lb.  Hence,  the  total  initiation 
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energy  in  the  tests  described  later  is  equal  to  the  blasting  cap  energy 
added  to  the  energy  of  the  Detasheet  used. 

4.  Experimental  Apparatus 

The  detonation  chamber  shown  in  Figure  30  is  mounted  on  a  Unistrut 
framework  with  the  chamber  centerline  approximately  4  feet  above  the 
floor.  The  fuel  storage,  motor -driven  fuel  pressure  pulse  generator, 
fuel  valving,  and  electrical  solenoids  are  also  mounted  within  the  Unistrut 
framework.  A  rack  cabinet  Immediately  to  one  side  of  the  detonation 
chamber  houses  most  of  the  associated  electrical  and  electronic  equipment. 
Figure  31  is  a  schematic  of  the  experimental  apparatus. 
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Figure  31.  Schematic  of  Experimental  Apparatus 


The  detonation  chamber,  without  the  breech-like  explosive  charge 
holder,  is  28  3/4  inches  along  the  centerline.  The  top  and  bottom  bars 
diverge  at  a  20°  total  Included  angle.  The  inside  dimensions  of  the 
chamber  are  2.05  inches  wide,  1  inch  high  at  the  narrow  end,  and  11 
inches  high  at  the  open  end.  The  top  and  bottom  bars,  as  shown  in 
Figure  32,  are  constructed  of  1  1/2  inches  x  4  inches  aluminum  stock  with 
1/2  inch  x  1  inch  grooves  milled  along  the  longitudinal  edges.  The  side 
windows  fit  into  these  grooves  and  are  held  in  place  by  larger  1/2-inch- 
thick  cold  rolled  steel  side  plates  bolted  through  the  aluminum  top  and 
bottom  bars  with  hardened  steel  high-strength  bolts.  For  the  initial  tests 
aluminum  plates  have  been  substituted  for  the  glass  windows.  However, 

1  inch-thick  Pyrex  windows  are  available  and  will  be  installed  when 
photographic  data  is  required. 

A  breech,  which  holds  the  detonator  cap  and  plastic  explosive,  is 
mounted  to  the  narrow  upstream  end  of  the  chamber.  Outward  details  of  the 
breech  are  apparent  in  Figure  33.  The  breech  is  constructed  of  3-lnch- 
dlameter  round  steel  stock  4  Inches  in  length.  In  the  initial  design,  the 
breech  Included  a  cylindrical  cavity  7/8  inch  in  diameter  and  2  1/2  inches 
deep  in  which  the  explosive  charge  was  placed.  One  end  of  this  cavity 
was  open  to  the  test  chamber;  at  the  bottom  (or  back)  of  this  cavity  a 
9/32-lnch-dlameter  hole  was  drilled  1/4 -Inch  deep  to  hold  the  end  of  the 
detonator  cap.  A  1/8 -inch -diameter  hole  was  bored  the  remaining  dis¬ 
tance  to  allow  the  detonator  lead  wires  to  be  fed  outside.  The  breech  was 
pressed  into  a  1 1/2  Inches  x  9  inches  square  aluminum  block  which  was 
bolted  to  the  chamber  backplate  and  sealed  wllh  an  O-rlng  seal. 

The  initial  blast  wave  <1.  e. ,  no  fuel)  tests  made  with  this  cylindrical 
cavity  in  the  breech  seemingly  indicated  that  the  design  was  satisfactory. 
However,  as  the  techniques  for  determining  wave  front  velocities  (using 


Figure  32.  Inside  of  Test  Chamber  Viewed  from  Exit  End 
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Figure  33.  Test  Chamber  with  Fuel  Sub-manifolds 
(top  and  breech  (left)) 


Figure  34.  Sub-manifold  with  Seven  0.008  in.  ID  Needles 
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multiple  pressure  switches,  etc.)  were  Improved  and  the  data  became  more 
reproducible,  an  anomaly  appeared.  The  Indicated  velocity  of  the  blast 
wave  showed  sudden  Increases  at  downstream  chamber  positions,  although 
blast  waves  In  a  diverging  chamber  would  normally  be  expected  to  slow 
down  with  Increasing  values  of  r. 

On  the  hypothesis  that  a  second  shock  was  overtaking  the  leading  shock 
(blast  wave)  and  thereby  reinforcing  and  Increasing  the  velocity  of  the  lead¬ 
ing  wave,  the  pressure  traces  (see  later  discussion  regarding  pressure 
transducers  usee)  were  examined  In  detail.  In  many  cases  It  was  apparent 
that  the  pressure  at  a  given  upstream  position  Increased  sharply  with  the 
leading  wave,  and  then  Increased  again  rather  sharply  a  few  microseconds 
later. 

Since  both  velocity  and  pressure  measurements  supported  the  thesis 
that  a  second  wave  (or  waves)  was,  presumably,  reflecting  from  the  backend 
of  the  cylindrical  explosion  cavity,  It  was  decided  that  this  cavity  should  be 
redesigned.  A  new  Insert  for  the  breech  assembly  was  therefore  made. 

The  explosive  charge  cavity  In  this  modified  breech  Is  a  compromise  be¬ 
tween  a  cylinder  and  a  cone  so  that  the  cavity  Is  flared  out  to  match  the 
1-lnch  by  2 -Inches  entrace  to  the  tee',  chamber  proper.  Also,  the  explo¬ 
sive  charge  can  now  be  located  very  nearly  at  what  would  be  the  apex  of 
the  top  and  bottom  sides  If  they  wer  ■  extended. 

Although  there  are  surely  sonr  reflected  waves  with  this  new  breech 
design,  the  problem  of  a  significant  second  wave  has  seemingly  disappeared. 
Essentially  all  of  the  blast/detonatlon  wave  test  results  reported  In  this 
report  were  obtained  with  this  Improved  breech. 

Fuel  Is  fed  Into  the  chamber  through  77  capillary-tube  needles.  The 
needles  are  In  11  groups  of  7  needles;  the  groups  are  2  5/32  Inches  apart 
and  the  needles  In  each  grotg>  are  1/4  Inch  apart.  The  stainless  steel 


needles  are  2  inches  long,  1/16  inch  outside  diameter,  and  .008  inch  inside 
diameter.  These  needles  are  soldered  into  sub-manifolds  of  3/16  inch 
stainless  steel,  seven  needles  to  each  sub-manifold  as  shown  in  Figure  34. 
The  sub-manifolds  are  connected  to  the  main  manifold  with  11  Conax  con¬ 
nectors.  The  main  manifold  is  constructed  of  a  1/2  inch  outside  diameter, 
1/4  inch  inside  diameter  brass  tube  which  runs  the  length  of  the  chamber 
top  bar.  Fuel  is  fed  into  the  lower  end  and  a  pressure  transducer  is 
connected  to  the  upper  end  through  a  solenoid  which  is  opened  during  the 
detonation  to  prevent  damage  to  the  transducer.  The  main  and  sub-mani¬ 
folds  were  constructed  so  that  gas  could  not  be^  trapped  in  any  connections 
or  fittings.  Continuous  gas  bleeds  are  provided  on  each  sub-manifold. 

Fuel  is  stored  in  a  one -cubic -foot  spherical  tank.  In  normal  operation, 
fuel  is  forced  out  by  pressurizing  the  tank  to  approximately  50  pslg  with 
nitrogen.  Fuel  flows  to  the  motor -driven  pressure  pulse  generator  through 
a  60-micron  filter,  a  solenoid  shut-off  valve,  a  rotometer,  and  a  throttling 
valve.  The  pulse  generator  Is  driven  by  a  1/3  horsepower  variable  speed 
motor.  Motor  speed  is  indicated  by  a  photo -transistor,  light  source,  and 
electronic  counter  arrangement  which  gives  1  pulse  output  for  48  pressure 
pulses  In  the  fuel. 

When  the  fuel  system  is  operated  in  the  setup  mode,  fuel  flows  into  the 
chamber  Into  a  Plexiglas  catch  pan,  and  a  pressure  transducer  is  con¬ 
nected  to  the  main  manifold.  After  the  proper  pressure  and  disturbance 
frequency  are  adjusted,  the  system  is  switched  to  the  ready  mode.  In  this 
mode,  fuel  is  switched  to  a  bypass  which  is  throttled  to  reflect  the  same 
pressure  and  mass  flow  as  the  needles  represent.  The  pressure  trans¬ 
ducer  is  also  connected  to  the  bypass  line.  In  the  automatic  firing  sequence, 
the  fuel  is  switched  from  the  bypass  to  the  main  manifold  2  1/2  seconds 
prior  to  detonation.  Hie  bypass  arrangement  causes  the  fuel  jets  and  thus 

the  droplets  from  the  needles  to  be  established  more  quickly. 

' 


To  prevent  fuel  from  wetting  the  floor  of  the  detonation  chamber  during 
the  2  1/2  seconds  the  drops  are  stabilizing,  a  1/4-inch-thick  absorbent  pad 
covers  the  entire  floor.  The  pad  is  snapped  out  approximately  80  msec  prior 
to  detonation.  In  this  way  a  minimum  amount  of  fuel  wets  the  chamber  floor. 
In  spite  of  these  precautions,  evidence  of  film  combustion  on  the  chamber 
walls  was  observed.  This  phenomenon  was  qualitatively  established  on  the 
basis  of  a  special  series  of  tests.  A  series  of  blast  wave  (no  fuel  adde<$ 
tests  were  performed  Immediately  following  a  series  of  detonation  tests 
without  cleaning  the  chamber  walls.  These  special  runs  revealed  both 
audible  and  visual  evidence  typical  of  kerosene  combustion.  Quantitative 
indications  of  film  combustion  were  also  obtained  from  detonation  tests  by 
a  comparison  of  the  data  taken  before  with  that  taken  after  careful  cleaning 
of  chamber  walls.  Such  comparisons  of  detonation  test  results  showed 
that  film  combustion  apparently  results  in  greater  scatter  in  the  experi¬ 
mental  data  as  well  as  a  modest  Increase  in  the  magnitude  of  wave  pres¬ 
sure  and  velocity.  These  Increases  were  more  evident  in  the  upstream 
portion  of  the  chamber  where  the  ratio  of  wall  surface  to  cross-sectional 
area  is  greatest. 

In  view  of  the  indication  that  residual  fuel  had  remained  on  the  chamber 
wall  and  reacted  in  later  tests,  a  better  procedure  for  chamber  cleaning 
between  tests  was  adopted.  The  chamber  walls  are  now  washed  out  with 
trichloroethylene  between  tests  involving  liquid  fuel.  The  effects  of  any 
reaction  with  a  film  of  residual  fuel  have  seemingly  been  eliminated. 

The  results  of  tests  involving  fuel  presented  in  this  report  were  obtained 
using  the  Improved  cleaning  procedure. 

The  test  results  presented  in  this  report  were  obtained  by  means  of 
three  pressure  transducers  and  several  pressure  switches.  The  output 
from  these  transducers  and  switches  were  recorded  using  three  Tektronix 
oscilloscopes,  each  with  an  oscilloscope  camera. 


One  pressure  switch  and  three  pressure  transducers  were  mounted  in 

the  bottom  of  the  chamber  0. 18,  0.  54,  1.  39,  and  2.  28  feet,  respectively, 

from  the  vertical  back  plate  of  the  test  chamber.  The  pressure  switch  is 

the  first  device  to  sense  the  blast  wave  from  the  explosive  charge.  The 

signal  from  this  switch  triggers  the  scopes  and  provides  a  point  of  time 

reference  for  all  other  data  recorded.  The  pressure  transducers,  shown 
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in  Figure  35,  are  similar  to  those  described  by  Rao  .  As  shown  in  this 
figure,  a  .  125>inch-diameter,  0. 0 50 -inch -thick  lead-zlrconate  titanate 
(PZT^  element  is  soldered  between  two  .  125 -inch -diameter  copper  rods, 

.  375  Inch  and  6.  50  inches  long.  Since  copper  has  approximately  the  same 
acoustic  impedance  as  the  PZT  crystal,  very  little  reflection  of  the  longi¬ 
tudinal  elastic  wave  occurs  at  the  Interfacial  surfaces.  The  elastic  wave 
does  reflect  off  the  far  end  of  the  copper  rod  and  affect  the  transducer  out¬ 
put  but  at  a  time  after  most  of  the  desired  pressure  data  has  been  recorded. 
PZT  exhibits  a  very  low  mechanical  Q,  and  this  reduces  problems  with 
oscillation  or  ringing  of  the  crystal.  The  rod  crystal  assembly  is  potted 
in  a  brass  housing  with  Dow  Corning  RTV  521  Silastic.  The  short  copper 
rod  Is  connected  electrically  to  the  outer  housing  by  a  very  fine  copper 
wire,  and  the  end  of  the  long  copper  rod  is  connected  to  the  center  of  a 
BNC  connector  attached  to  the  end  of  the  transducer  housing.  A  typical 
calibration  curve  for  one  of  the  transducers  Is  shown  In  Figure  36.  These 
transducers  were  individually  calibrated  in  a  detonation  tube  using  known 
gaseous  detonation  conditions. 

In  addition  to  the  three  pressure  transducers  and  the  one  pressure 
switch  referred  to  above,  eight  pressure  switches  have  been  installed  In  a 
radial  line  along  one  of  the  side  walls  of  the  test  chamber.  The  distance 
between  each  pressure  switch  Is  2. 4  Inches.  These  switches  are  connected 
Into  a  raster  circuit  so  that  as  the  leading  wave,  blast  or  detonation,  passes 
each  switch  in  succession,  the  event  Is  Indicated  on  a  scope  trace. 


L  Copper  a  125  inch 
Diameter  by  0.375  inches 

2.  Lead  Zirconate  Tftanate 
(PZT)  0.125  Inch 
Diameter  by  (L  050  inches 
Thick 

3.  Copper  Ql  125  inch 
Diameter  by  6. 75  inches 
Long 

4  Silicone  Rubber 

5.  Brass 


Flcurt  M.  Pleasure  Tranaduar  Asaemtaiy 


Figure  36.  Pressure  Transducer  Calibration 


These  pressure  switches  are  made  up  essentially  of  a  1/3  inch  stain¬ 
less  steel  tube  with  an  insulated  center  conductor  and  a  thin  diaphragm 
soldered  across  the  end  of  the  tube  exposed  to  the  chamber.  Figure  37 

is  a  sketch  of  one  of  these  switches.  Although  pressure  switches  of  a 
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basically  similar  nature  have  been  used  at  this  laboratory  for  years  in 
detonation  work,  the  switches  used  here  are  of  a  new  design.  The  main 
virtues  of  this  revised  design  over  switches  previously  used  are  its 
small  diameter  and  extended  length.  Since  a  large  number  of  switches 
were  to  be  mounted  fairly  close  together  In  a  fairly  thick  wall,  this  par¬ 
ticular  design  appeared  to  be  appropriate. 

Prior  to  the  design  and  fabrication  of  these  pressure  switches  several 
small  and  Inexpensive  Pin-ducers*  had  bejn  Installed  in  the  same  locations. 
While  these  inexpensive  pin-ducers  apparently  had  the  desired  character¬ 
istics  of  fast  response  time  and  relatively  high  outputs  (for  the  pressure 
increase  expected),  they  failed  to  survive  a  large  number  of  blasts.  After 
a  few  tests  they  generally  became  erratic,  apparently  due  to  erosion  of 
the  very  thin  coating  on  the  pressure -sensitive  end.  When  it  was  learned 
that  these  Pin-ducers  were  really  designed  to  survive  only  one  strong  blast 
wave,  no  further  attempts  were  made  to  use  them.  Effort  was  therefore 
directed  toward  the  development  of  the  pressure  switches  discussed  above. 

5.  Experimental  Procedure 

A  systematic  experimental  procedure  was  established  to  minimize  the 
chances  of  wasted  runs,  either  as  the  result  of  equipment  failure  or  human 
error.  The  sequence  of  events  prior  to  a  detonation  run  is  as  follows: 


♦These  Pin-ducers  are  co-axial  piezoelectric  (modified  barium  tltanate) 
transducer  pins,  Channel  Industries  part  number  PK 14-X. 
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Copper,  0.0625  inch 
diameter  rod, polished  end. 

Stainless  Steel,  0. 125  inch 
OD  tubing. 

1/2  mil,  Stainless  Steel 
Diaphragm,  soldered  to 
tubing. 

Bakelite 


1.  Inspect  detonation  chamber  and  remainder  of  basic  hardware. 

2.  Inspect  and  adjust  pressure  and  time  of  arrival  measurement 
instrumentation. 

3.  Clean  detonation  chamber  thoroughly. 

4.  Install  fuel  catch-pan  In  chamber  and  start  fuel  flow. 

5.  Adjust  tael  pressure  and  pressure  pulse  generator  frequency  to 
operating  conditions. 

6.  Switch  fuel  system  to  ready  mode  and  remove  catch-pan. 

7.  Install  absorbent  pad  and  arm  pad  retraction  mechanism. 

8.  Install  detonator  cap  with  appropriate  amount  of  plastic  explosive 
attached  in  breech. 

9.  Secure  breech  to  backplate  studs. 

10.  Arm  detonator  circuit. 

11.  Arm  automatic  run  sequencer. 

12.  Open  scope  camera  shutters. 

The  start -run  switch  activates  warning  lights  and  an  eight -second 
duration  warning  horn  in  the  test  area.  Ten  seconds  after  start,  a  mechan 
ical  microswitch  timer  is  started.  The  timer  turns  fuel  on  immediately, 

2  1/2  seconds  prior  to  detonation  and  turns  off  within  1  second  after  detona 
tion.  Approximately  80  msec  prior  to  detonation,  the  absorbent  pad  is 
snapped  from  the  chamber. 

Immediately  following  detonation,  fuel  is  shut  off,  the  detonator  and 
automatic  run  sequencer  arc  disarmed,  and  the  photographs  of  the  oscillo¬ 
scope  traces  are  developed.  When  detonation  in  an  oxygen-enriched 
atmosphere  is  desired,  the  chamber  is  first  purged  of  air  by  filling  the 
chamber  with  premixed  oxygen  and  nitrogen  through  a  port  in  the  bottom 
near  the  pressure  switch  and  a  thin  film  of  plastic  is  taped  across  the 
open  exit  end.  Throughout  the  run,  the  premixed  N„-00  is  blown  slowly 
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into  the  chamber  to  maintain  the  desired  chamber  atmosphere.  The  pro¬ 
cedure  required  to  carry  out  blast  wave  runs  eliminates  steps  4  through  7. 


C.  RESEARCH  RESULTS  AND  DISCUSSION 
1.  Introduction 

It  is  unlikely  that  a  complete  and  detailed  analytic  solution  is  possible 
to  the  complex  problem  of  fuel-air  explosions.  The  formulation  and  solu¬ 
tion  of  such  a  complex  model  must  then  rely  heavily  upon  experiment. 
Theoretical  consideration  of  this  problem  was  initially  directed  toward  a 
detailed  discussion  of  the  two  limiting  solutions  of  a  complete  FAE  model — 
the  cylindrical  blast  wave  and  cylindrical  Chapman -Jouguet  detonation  wave. 

Experimental  FAE  research  has  similarly  been  directed  toward  the 
study  of  these  two  limiting  cases.  To  this  juncture,  work  has  focused 
upon  the  evaluation  of  detonation  initiation  and  propagation  characteristics 
in  a  fuel  cloud.  To  properly  achieve  this  goal,  thorough  experimentation 
was  required  for  both  the  blast  wave  as  well  as  the  detonation  wave.  Con¬ 
siderable  emphasis  has  been  placed  on  the  attainment  of  good  blast  wave 
data.  The  reasons  for  this  are  basically  three -fold.  First,  it  was  an 
expedient  way  to  debug  much  of  the  apparatus  and  Instrumentation  and  to 
gain  experience  and  confidence  in  the  operation.  Second,  It  was  deemed 
important  to  assess  the  degree  to  which  cylindrical  blast  waves  were  being 
generated.  Third,  good  documentation  and  understanding  of  the  blast  wave 
data  should  provide  an  experimental  technique  for  the  determination  of 
initiation  requirements  for  the  two-phase  detonation. 

Table  n  gives  an  abbreviated  summary  of  experimental  tests  performed 
during  the  course  of  this  research.  Work  progressed  from  the  initial  stage 
of  hardware  development  and  calibration  through  the  preliminary  testing 
phase  to  the  final  data  seeking  phase. 

The  information  obtained  during  the  controlled  tests  consisted  of  the 
wave  time  of  arrival  at  eight  different  chamber  radii  and  wave  pressure 
at  three  different  radial  positions.  The  wave  time  of  arrival  data  is 


TABLE  n.  SUMMARY  OF  EXPERIMENTAL  TESTS  PERFORMED 


Run  No. 

Type 

Explosive 

Charge1 

Comments 

• 

Hardware 

Development 

and 

Calibration 

Items  tested  included  pin 
transducers,  pressure 
switches,  pressure  trans¬ 
ducers  and  associru  jd 
instrumentation. 

1-30 

Blast  Wave 

1100-8400 

Original  breech  design 

31-41 

Blast  Waveb 

1100-8400 

Improved  breech 

42-46 

Detonation 

1100-7000 

O2  enriched  atmosphere 

47-60 

Detonation 

2500-8400 

Air  detonation  with  wall 
wetting  problem 

61-63 

Special 

2500 

Film  combustion  check 

64-76 

Blast  Waveb 

2500-8400 

Improved  breech 

77-85 

Detonatlonb 

2500-8400 

Air  detonation  with 
minimized  wall  wetting 

aTotal  energy  release  In  ft-lbf  as  computed  by  technique  of  Section  B.  3. 
Results  of  these  tests  presented  in  detail  In  the  discussion. 
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obtained  from  pressure  switches  positioned  in  the  chamber  side  wall. 

The  output  from  these  switches  is  displayed  by  means  of  a  raster  circuit 
on  an  oscilloscope.  An  example  of  time  of  arrival  data  is  given  in  Figure 
38  in  the  form  of  a  scope  picture  of  a  raster  trace.  The  wave  speed  of 
propagation  can  be  obtained  approximately  by  knowing  the  distance  between 
pressure  switches  and  the  time  between  raster  spikes.  However,  it  is 
determined  more  accurately  from  the  slope  of  the  smoothed  position-time 
data.  The  majority  of  the  data  discussed  in  the  following  sections  is 
presented  in  terms  of  a  time  measured  from  the  scope  trigger  time,  t\ 
Figure  39  is  representative  of  wave  pressure  data  taken  during  the  tests. 


The  following  sections  present  and  discuss  the  results  of  good  experi¬ 
mental  blast  wave  and  detonation  tests.  The  problem  of  relating  the  wave 
propagation  characteristics,  measured  experimentally  with  respect  to 
time  t',  to  a  time  associated  with  the  hypothetical  origin  of  the  cylindrical 
model  is  recognized  and  discussed.  Experimental  blast  wave  data  Is 
presented  and  briefly  discussed  in  terms  of  its  fit  to  the  strong  blast  wave 
theory. 


The  Importance  of  the  critical  blast  wave  radius  is  discussed  in  light 
of  both  experimental  and  analytical  detonation  work.  Experimental  detona¬ 
tion  and  blast  wave  results  are  compared  to  establish  the  probable  occur¬ 
rence  of  a  fuel-air  detonation. 

2.  Blast  Wave  Results 

Strong  shock  waves  (blast  waves)  were  produced  with  various  ignition 
energy  levels  and  without  any  fuel  present  in  the  chamber.  As  mentioned 
above,  the  data  were  obtained  by  measuring  the  time  of  arrival  of  the 
blast  waves  at  various  stations  along  the  chamber.  This  time,  t\  was 
measured  from  the  farthest  upstream  pressure  switch  which  was  used  as 
a  trigger  for  the  oscilloscopes.  Comparisons  of  experimental  data  obtained 
for  three  different  energy  levels  are  shown  in  Figures  40  and  41.  The 
energy  level  indicated  represents  the  grams  of  Deta sheet  used  in  addition 
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Figure  38.  Detonation  Run  Raster  Trace 


Figure  39.  Blast  Wave  Run  Pressure  Trace 
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Figure  40.  Blast  Wave  Data  Plotted  as  r  Versus  t’ 


Figure  41.  Blast  Wave  Data  Plotted  as  M  Versus  r 


to  the  blasting  cap.  Four  runs  were  made  at  each  energy  level.  Figure  40 
shows  the  shock  wave  position  versus  t',  and  the  curve  drawn  for  each 
energy  level  is  the  mean  curve  for  the  4  runs.  Figure  41  shows  the  same 
data  plotted  as  shock  Mach  number,  M,  versus  r.  Values  of  M  were 
derived  from  M  =  v/a  where  v  is  shock  velocity,  a  is  the  speed  of  sound, 
and  v  was  obtained  from  the  slope  of  the  r  - 1*  curves.  As  seen,  the  data 
are  quite  reproducible,  and  the  scatter  is  reasonably  small.  As  expected, 
the  lower  energy  levels  produced  weaker  blast  waves.  Near  the  end  of  the 
chamber  the  shock  strengths  have  decreased  below  M  =  4  for  each  energy 
level  and  hence  would  begin  to  deviate  from  strong  blast  wave  behavior. 

In  all  cases,  the  shocks  gradually  decay  towards  Mach  waves. 

It  was  of  interest  to  Interpret  these  results  in  light  of  cylindrical  blast 
wave  theory.  In  order  to  do  this,  it  was  necessary  to  determine  both  the 
theoretical  geometric  origin  and  the  time  origin  of  the  modeled  line-source 
explosion.  First  of  all,  a  radius  value  of  r  =  0  was  assumed  to  be  at  the 
hypothetical  apex  of  the  chamber;  that  is,  the  assumed  origin  of  the  cylin¬ 
drical  wave  the  equipment  was  designed  to  model.  This  is  shown  in 
Figure  42. 

I 
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Figure  id.  Origin  of  FAB  Model 


Once  this  geometric  origin  was  established,  it  was  necessary  to  deter¬ 
mine  the  time  it  would  take  a  blast  wave  to  travel  from  the  origin  to  the 
trigger  station.  This  time  is  represented  by  tQ.  As  all  of  the  time  values 
had  previously  been  referenced  from  the  trigger,  the  time  value,  t,  refer¬ 
enced  from  the  origin,  is  determined  by 

t  =  t’-tQ  (138) 

where  t  was  also  measured  from  the  trigger  (i.e. ,  t  was  negative).  This 
is  also  illustrated  in  Figure  42.  Using  the  experimental  data  graphs  of  r 
versus  t'  for  blast  waves,  a  value  of  t  was  determined  using  a  geometrical 
technique  described  by  T. R.  Running^®).  The  value  of  t  arrived  at  for  an 
initiation  energy  level  produced  by  2. 5  gm  of  Detasheet  "C"  and  a  blasting 
cap  was  -45  psec.  This  is  fairly  close  to  the  value  of  -40  jisec  obtained  by 
extrapolating  the  t'  versus  r  curve  to  the  value  of  r  -  0. 

With  the  origin  determined,  it  was  desired  to  establish  an  equation 
that  would  describe  the  experimental  results.  The  form  of  this  equation 
was  naturally  chosen  to  be  that  of  a  strong  blast  wave,  that  is: 

t  =  t'  -  t  =  b  rc  (139) 

i 

where  b  and  c  are  constants.  Taking  the  logrithm  of  both  sides  of  this 
equation, 

log  (t*  -  tQ)  =  log  b  ♦  c  log  r  .  (140) 

This  is  the  equation  of  a  straight  line  in  slope -Intercept  form  in  which  log  b 
is  the  intercept  on  the  log  (t*  -  tQ)  axis  and  c  is  the  reciprocal  of  the  slope 
when  log  r  is  plotted  against  log  (r  -  tQ).  However,  upon  plotting  the 
experimental  results,  a  straight  line  was  not  obtained,  although  a  smooth 
curve  did  result.  This  curvature  could  be  altered  and  eliminated  by  employ 
Ing  different  values  for  tQ.  Thus,  by  using  a  series  of  successive  approxi¬ 
mations,  a  new  value  of  t  was  determined  which  produced  a  straight  line 


in  the  log  (t’  -  t  )  versus  log  r  graph.  This  is  illustrated  in  Figure  43 

which  shows  plots  obtained  using  different  values  of  tQ.  The  sensitivity 

of  the  resultant  curve  to  the  value  of  t  should  be  noted.  The  value  of  t 

o  o 

derived  for  the  straight  line  was  -25  jisec. 


The  value  of  c  was  then  determined  as  the  reciprocal  of  the  slope  of 
the  above  mentioned  graph, 


A  log  (t*  -  tQ) 
c  “  A  log  r 

The  value  of  c  for  the  energy  level  discussed  was  found  to  be  1. 56. 


(141) 


A  value  of  b  was  ascertained  in  two  ways,  first  by  finding  the  line- 
intercept  value  [on  the  log  (t*  -  tQ)  axis]  of  the  straight  line  graph,  and 
second,  by  picking  a  point  on  the  graph,  using  the  values  of  c  and  t 
already  determined  and  then  using  the  empirical  equation 


t  =  t'  -  t  =  brc 
o 


(142) 


The  former  method  yielded  a  value  of  b  =  2  and  the  latter  b  =  2.03. 


Using  these  values  for  the  constants,  the  derived  equation  for  the 
data  (2. 5  grams)  was  found  to  be: 


t  *  t'  +  25  »  2r1- 


(145) 


where  t  and  t'  were  measured  in  microseconds  and  r  in  inches.  This  equa 
tlon  Is  plotted  in  Figure  44,  along  with  the  experimental  data  plot  in  which 
a  value  of  tQ  «  >45  n sec  was  used  (the  value  originally  obtained  from  the 
experimental  data). 


According  to  strong  Mast  wave  theory,  the  relationship  between  wave 
radius  and  time  is: 


t .  (p/K)1/2  r2 


(144) 
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Figure  44.  Plot  of  r  Versus  t  for  Experimental  Values  and  Derived  Blast  Wave  Equation 


It  is  lo  be  noted  that  the  experimentally  derived  value  of  e  is  1.  56  as  com¬ 
pared  to  the  pure  theoretical  value  of  2.0.  Further,  assuming  no  energy 

1/2 

losses  (such  as  losses  to  the  chamber  walls),  a  rough  value  of  (p/E)  ' 
was  calculated  to  be  0.  338.  This  value  corresponds  to  the  value  of  b  =  2 
in  the  experimentally  derived  equation.  Some  discrepancy  between  the 
strong  blast  wave  theory  and  the  experimental  equation  could  have  resulted 
from  the  fact  that  near  the  end  of  the  chamber  a  strong  blast  wave  s  uo 
longer  experienced  (l.  e. ,  one  in  which  the  pressure  ratio  is  much  greater 
than  1.0).  Other  explanations  for  the  discrepancy  include  the  fact  that  the 
charge  is  not  a  true  line  source  of  energy  addition  and  also  that  the  energy 
release  is  not  Instantaneous.  The  complete  Interpretation  of  the  blast 
wave  data  is  not  yet  complete,  and  these  points  of  discrepancy  will  be 
considered  further. 

3.  Detonation  Test  Results 

Experimental  detonation  runs  were  performed  at  three  different  initia¬ 
tion  energy  levels  with  kerosene,  C  drops  (approximately  360  microns) 
uniformly  distributed  throughout  the  chamber.  The  overall  fuel -air  mixture 
ratio  for  these  tests  was  essentially  stoichiometric,  the  equivalence  ratio, 
ft,  was  1.04. 

The  results  of  these  tests  are  presented  as  raw  data  plots  of  chamber 
radial  distance,  r,  versus  time,  t’,  for  the  three  energy  levels  examined. 
Again  a  number  of  runs  were  made  at  each  energy  level,  and  the  curves 
drawn  are  the  mean  of  the  Individual  runs.  The  data  is  further  reduced  by 
taking  dr/dt  of  these  curves  to  derive  wave  Mach  number,  M,  and  plotting 
this  against  chamber  radius.  Figures  45,  46,  and  47  are  plots  of  r  versus 
t  for  initiation  energy  levels  of  8370  (2. 5  grams),  5470  (1. 5  grams),  and 
2561  (0. 5  grams)  ft-lbf,  respectively.  These  energies  are  computed 
directly  from  a  known  quantity  of  Datasheet  and  represent  total  energy 
releases  of  the  Initiating  charge,  including  the  blasting  caps. 


t'  -  Time  Measured  from  Trigger,  jj  sec. 

Figure  46.  Detonation  Data  Plotted  as  r  Versus  t’  for  1.  5  gm  Detasheet 


Figure  47.  Detonation  Data  Plotted  as  r  Versus  t'  for  0.  5  gm  Detasheet 


Figure  45  reveals  very  strong  agreement  between  four  successive  high 
energy  experimental  runs,  a  very  welcome  and  rarely  achieved  experimental 
result.  One  can  clearly  see  that  downstream  of  the  mid-chamber  point  the 
data  takes  on  a  linear  relationship  between  r  and  t\  This  result  is  of 
extreme  interest  since  it  is  this  precise  wave  behavior  which  is  commonly 
associated  with  propagation  of  a  sustained  Chapman -Joqguet  detonation. 
Examination  of  Figures  45  and  46  reveals  the  same  phenomenon  oeoin  mg 
but  with  slightly  varying  properties.  The  curve  of  Figure  46  also  .as 
small  data  scatter  and,  significantly,  the  linear  portion  of  the  curvo  is 
clearly  longer  than  in  Figure  45.  This  latter  fact  is  of  interest  when  the 
subject  of  critical  blast  wave  radius  is  considered.  Figure  47  reveals  that, 
for  lower  energies  than  those  used  in  runs  for  Figures  45  and  46,  the  data 
scatter  is  greater  and  there  is  evidence  that  the  wave  may  be  slowing  down. 

Figure  48  is  a  plot  of  Mach  number  versus  chamber  radius  for  the 
same  three  initiation  energies.  The  results  are  quite  graphic  and  descrip¬ 
tive  of  the  process  occurring.  Before  discussing  Figure  48  further,  calcu¬ 
lations  were  performed  using  Equation  (1)  to  seek  the  critical  blast  wave 
radius  as  a  function  of  initiation  energy  for  our  particular  simplified  FAE 
model.  For  cylindrical  geometry  Equation  (1)  becomes 


where  all  the  definitions  remain  unchanged.  The  combustion  energy  release 
in  the  detonation  wave  is  needed  to  make  this  calculation.  Now  it  is  known 
that  use  of  the  constant  pressure  heat  of  combustion  would  lead  to  consider¬ 
able  Inaccuracy.  Consequently,  it  is  better  to  take  the  experimentally  deter¬ 
mined  Mach  number  of  detonation  (or  more  precisely,  the  pressure  immedi¬ 
ately  behind  the  shock)  and  then  calculate  Q  on  the  basis  of  simplified  one- 
dimensional  Chapman -Jouguet  detonation  theory.  Within  these  limitations 
the  conservative  equations  of  energy,  mass,  and  momentum  may  then  be 
written 


Figure  48.  Detonation  Plotted  as  M  Versus  r 
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hi+“T=h2+T‘"Q 


P1V1  =  P3V3 


(145) 


where 


(3) 


P1  +  plvl2  =  p3  +  p3v32 


(2) 


C-J 


Shock 


Mt  (1) 


For  a  calorlcally  perfect  gas, 


y  - 1  p 


(146) 


Equation  (145)  may  be  rewritten  and  solved  for  Q  in  terms  of  , 

Yv  and  Pj,  P3,  y3;  thus 


0  1  -  ./I  1\  y3  p3  'l  **1 

^  *  2  P1  "  P3^p3  +  Pi/  Vg  +  1  Pg  ”  -  1  Pj 


1  pl 


(147) 


Knowing  that 


Pt  »  2116.  2  lbf/ft' 
Tt  -  537  °R 
rx-  1.4 
y3  =  1.2 


(148) 


the  remainder  of  the  parameters  must  be  known  in  order  to  solve  for  a 
numerical  value  for  Q.  It  is  known  that  for  nearly  all  detonations 


1.80 


and 


p3 

— -  z  0.  55 
P2 

To  find  (Pg/  Pj)  we  have 

pS  p3  p2  .  „  p2 
pl  p2  pl  pl 


(149) 


(150) 


At  the  pressure  transducer  position  (M),  the  experimental  shock  pressure 
ratio  Pg/Pj  1128  recorded  as 


(151) 


for  both  the  2. 5-g  'am  and  l.  5-gram  energy  runs.  The  value  of  Equation 
(151)for  ( Pg/ 1  )  was  assumed  for  calculations  and  Is  justified  as  It  Is  the 
ratio  In  the  region  nearest  the  point  the  wave  begins  to  propagate  as  a 
C-J  detonation.  Thus,  on  this  basis 

—  =  12. 55  .  (152) 

Pl 

The  ambient  density,  was  computed  from  the  perfect  gas  law. 

The  combustion  energy  release  of  Equation  (147)  was  found,  upon  a 
substitution  of  the  above  values,  to  be 

Q«  8.975  x  105  ft-lK/lbm  .  (153) 

A  value  of  critical  energy  release  may  then  be  computed  by  substitution  of 
the  values  of  Q,  p,  and  the  desired  energy  level  into  Equation  (144).  The 


desired  energy  level  in  ft-lbf/ft  is  arrived  at  by  taking  a  figure  for  total 
initiation  energy  release,  computed  by  the  technique  of  Section  B.  3,  and 
correcting  it  by  a  factor  established  by  the  geometry  of  the  selected  FAE 
model. 

The  calculated  energy  release  must  be  related  to  the  model  geometry 
by  considering  that  it  was  released  within  a  20°  segment  and  over  a  length 
of  2.05  inches  of  the  modeled  cylinder.  Therefore,  the  effective  em  gy 
release,  Eq,  required  to  compute  r^  with  the  calculated  energy  release 
and  the  model  geometry  can  be  related  by 


Figure  49  displays  the  results  of  plotting  critical  blast  wave  radius  as  a 
function  of  initiation  energy  for  the  chosen  FAE  model. 

It  is  now  of  interest  to  return  to  Figures  45  and  48  and  pick -off  experi 
mental  values  of  r,  for  comparison  with  the  computed  values.  The 
results  are  given  in  Table  m  where  the  total  initiation  energy  (ft-lbf)  has 
been  multiplied  by  the  geometric  correction  factor  to  yield  Eq  in  ft-lbf/ft. 


TABLE  m.  COMPARISON  OF  r*..  AND  r* 

*th  *ex 


Explosive 

Charge 

-  ft-ltaf 
o  ft 

r*th 

r*ex 

Per  cent 
Difference 

2. 5  gm 

8. 82  x  105 

25.6 

20.0 

21.8 

1. 5  gm 

5.77xl05 

20.7 

18.0 

13.0 

0. 5  gm 

2.70xl06 

14.2 

22.0 

54.8 

The  results  appear  very  favorable,  yielding  an  average  difference  of  only 
17  per  cent,  If  the  lowest  energy  level  could  be  disregarded. 


A  comparison  of  experimental  blast  wave  and  detonation  results  is 
provided  by  Figures  50  and  51.  The  differences  on  an  r-t  plot  are  quite 
subtle.  However,  it  is  the  slope  of  the  r-t  curve  which  is  of  interest,  and 
it  can  be  noted  that  the  slope  of  the  blast  wave  curve  continues  to  decrease, 
indicating  a  decreasing  wave  velocity.  This  is  not  the  case  on  the  uetona- 
tion  curve  beyond  a  certain  point.  On  an  M-r  plot  the  differences  are 
much  more  graphic.  Both  the  blast  wave  and  the  blast  wave  initiated 
detonation  begin  at  closely  similar  Mach  numbers,  as  they  should.  As 
the  waves  propagate  with  increasing  radius,  the  blast  wave  decays  markedly. 

It  is  clearly  evident  from  the  foregoing  that  a  wave  structure  was 
propagated  with  characteristics  normally  associated  with  those  of  a 
Chapin  an -Jouguet  detonation.  While  the  results  appear  most  promising, 
further  stu^  l®  required  to  establish  greater  experimental  documentation 
of  the  propagation  of  cylindrical  detonation  waves.  The  experimental 
'  results  must  also  be  further  compared  with  theoretical  predictions. 


Transducer 


Figure  50.  Comparison  of  Experimental  Detonation  and  Blast  Wave  Data 

on  an  If  Versus  r  Plot 


Figure  51.  Comparison  of  Experimental  and  Blast  Wave  Data 

on  an  M  Versus  r  Plot 


SECTION  IV 
FUTURE  PLANS 

A.  THEORETICAL  CONSIDERATIONS 

During  the  coming  year  It  is  planned  to  continue  our  theoretical 
analysis  of  FAE  as  outlined  below. 

1.  C-J  Waves  with  Cylindrical  and  Spherical  Symmetry 

Plans  Include  the  completion  of  the  computer  program  for  computing 
the  propagation  of  cylindrical  and  spherical  self-annular  C-J  waves.  This 
program  will  be  used  In  conjunction  with  the  analysis  described  below. 

2.  The  Influence  of  Droplets  on  Propagation 

As  indicated  in  this  report,  the  presence  of  the  fuel  in  the  form  of 
liquid  drops  influences  the  propagation  velocity  and  the  change  in  proper¬ 
ties  across  a  C-J  wave.  It  is  planned  to  incorporate  the  revised  C-J 
conditions,  both  in  the  self  similar  C-J  calculations  and  in  the  NASA 
program,  for  computation  of  C-J  velocities  in  hydrocarbon  fuel.  Efforts 
will  also  be  made  to  develop  simple  analytical  expressions  which  show 
the  Influence  of  the  droplets  on  the  C-J  conditions. 

3.  Idealized  FAE  Calculations 

The  Idealized  FAE  consists  of  a  blast  wave  to  r  *  r*,  followed  by  a 
C-J  wave  for  r  >  r„.  This  idealized  model  shall  be  used  to  study  the 
behavior  of  FAE  for  various  blast  energies  and  for  different  fuels  and 
equivalence  ratios. 

4.  Ground  and  Dynamic  Impulse 

The  Impulse  calculations  shall  be  extended  as  described  in  the  present 
report  with  the  objective  of  comparing  the  effectiveness  of  the  FAE  to 
blast  waves.  Some  of  the  constants  appearing  in  the  theory  presented 
here  also  remain  to  be  calculated.  The  Influence  of  shock  diffraction 


upon  dynamic  impulse  will  also  be  examined.  A  number  of  investigators 
have  studied  the  Interaction  of  shock  waves  with  various  objects,  and  it  is 
planned  to  apply  these  results  to  the  FAE. 

5.  The  Effect  of  Side  Relief 

Side  relief  due  to  the  finite  size  of  the  fuel  cloud  will  have  an  important 
effect  on  the  pressure  and  velocity  distribution  within  a  FAE.  An  attempt 
was  made  to  develop  simple  models  which  will  Indicate  the  major  effects 
of  finite  cloud  size. 


6.  Initiation 

The  problem  of  blast  initiation  of  FAE  will  be  studied  further.  In 

particular,  an  attempt  will  be  made  to  compare  the  expression  for  (E  ) 

°  crit 

presented  in  this  report  with  experimental  measurements.  The  effect  of 
cloud  geometry  upon  the  initiation  process  will  also  be  Investigates 


7.  Non  Uniform  Droplet  Distribution 

A  similarity  solution  for  self  similar  C-J  waves  with  a  non-uniform 
distribution  of  fuel-oxidizer  ratio  has  already  been  presented.  This  and 
other  similarity  solutions  will  be  examined  to  assess  the  main  effects  of 
the  non-uniform  distribution  of  droplets  within  the  fuel  cloud. 


B.  EXPERIMENTAL  RESEARCH 


The  experimental  facility  and  Instrumentation  appear  to  be  working 
very  satisfactorily.  A  fair  amount  of  blast  wave  data  (no  fuel)  has  been 
obtained,  and  the  results  are  quite  reproducible  with  relatively  small 
data  scatter.  Interpretation  of  the  data  in  terms  of  theoretical  cylindrical 
blast  wave  behavior  presents  some  difficulty.  Accordingly,  it  is  planned 
to  further  reduce  and  analyse  the  present  data.  If  the  results  indicate  the 
need,  some  further  runs  would  be  made. 
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The  experiments  conducted  in  kerosene -air  have  been  extremely 
encouraging.  With  sufficiently  high  initiation  energies,  the  blast  wave 
decays  in  strength  until  the  critical  radius  where  the  detonation  Mach 
number  is  realized.  Evidently,  beyond  this  point  a  steady  state  Chapman - 
Jouguet  detonation  is  established.  At  lower  energy  levels  there  is  evi¬ 
dence  that  the  wave  is  not  self-supporting  and  hence  would  gradually 
quench.  This  facet  is  of  considerable  importance  to  the  initiation  require¬ 
ments,  and  it  is  planned  on  pursuing  this  point  in  some  detail  in  order  to 
establish  a  criterion  for  initiation.  It  is  planned  that  more  experimental 
data  and  analysis  of  the  data  will  be  obtained  in  the  Immediate  future. 

Upon  achieving  satisfactory  results  on  the  above,  detonation  experi¬ 
ments  would  then  be  conducted  for  a  few  different  kerosene -air  mixture 
ratios.  The  Influence  of  equivalence  ratio  on  initiation  energy  would  be 
noted.  After  this,  another  fuel,  with  properties  substantially  different 
from  those  of  kerosene,  would  be  tested.  Other  work,  such  as  various 
drop  sizes,  non-uniform  fuel  distribution,  and  high  speed  photography, 
would  be  Incorporated  upon  achieving  the  desirable  level  of  understand¬ 
ing  of  the  foregoing.  It  is  also  planned  to  add  a  receiver  section  to  the 
)* 

outlet  of  the  test  chamber  in  order  to  reduce  the  severe  noise  which 
accompanies  the  detonation  tests. 
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The  theoretical  aspects  of  cylindrical  strong  blast  waves  and  Chapman - 
Jouguet  (C-J)  detonations  are  treated  in  detail.  A  critical  radius  is  discussed 
which  divides  the  cloud  into  an  inner  blast  wave  zone  and  an  outer  detonation  zonql 
Expressions  for  the  evaluation  of  ground  and  dynaaic  iapulse  are  presented,  and 
the  propagation  of  a  C-J  detonation  through  a  cloud  of  fuel  drops,  the  initiation 
problem,  and  deviations  from  the  ideal  cylindrical  blast  aodel  are  considered. 

A  modified  computer  program  for  calculating  detonation  velocity  of  complex  hydro¬ 
carbon  fuels  is  discussed.  The  design  and  operation  of  an  experimental  facility 
to  study  the  propagation  of  two-phase  detonations  in  a  cylindrical  segment  of  a 
cloud  is  described.  Experiments  conducted  with  kerosene  drops  in  air  showed  that 
at  saall  radius  the  cylindrical  wave  decayed  as  a  shock  wave,  but  beyond  a  critica 
radius  the  wave  becomes  a  constant  velocity  two-phase  Chapman -Jouguet  detonation. 
The  experimentally  determined  critical  radius  agreed  quite  well  with  the  theoreti¬ 
cal  predictions.  The  results  lend  encouragement  to  the  prediction  of  threshold 
energy  levels  required  for  detonation  initiation  in  various  geometries. 
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